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SYNOPSIS 

SOME STUDIES ON STATISTICAL ASPECTS OE SIZE EFFECTS 
ON STRENGTH AND FRACTURE BEHAVIOUR OF MATERIALS AND 
FRACTURE RESISTANT DESIGN 

A Thesis Submitted 

In Partial fulfilment of the Requirements 
for the Degree of 

DOCTOR OF PKELOSOPHY 

"by 

C.V.S. ICAMSS VARA RAO 
Indian Institute of Technology 
Kanpur- 1 6 , India 


The classical approach of continuum mechanics 
of materials is inadequate to account for the phenomena 
like size effects on strength, ductile-brittle fracture 
transition etc. This is because the approach does not 
consider the existence of inherent structural defects 
like flaws in materials. Applying the theories of 
fracture mechanics and composite materials, the concepts 
of classical mechanics of materials are generalized 
to interpret size effects on strength, stiffness and 
ductile-brittle transition. Also, revised design 
methodologies' are evolved to ensure certain levels of 
safety from brittle fracture, with suitable optimality 
criteria. The present thesis is a contribution to the 
aforementioned areas. 


A critical review of the studies on size effects 



on strength of various materials and statistical theories 
of strength "based on weakest link concept and their 
application to various materials, is made. Following 
the observations in the review, it is emphasized that 
fitting test date, suitable to an a prior i chosen form 
of strength distribution function, as is done in present 
approaches, is limited in its application to a restricted 
class of material. s and is at best an approximation. 
Therefojre, a more general approach is proposed in this 
work by which a. strength distribution function can be 
obtained corf esc onding to an actual form of size-mean 
strength relation (obtained experimentally) of a 
particular material. Using this formulation and method 
of solution, it is shown that ^eibull’s distribution 
function can be obtained as a particular case corres- 
ponding to a particular form of size-mean strength 
relation. Construction of strength distribution 
functions corresponding to more general forms of size- 
mean strength relations needs resort to numerical 
methods which are outlined in this investigation. The 
approach developed herein to study size effects on 
strength of flaw sensitive materials is motivated by 
the work of Tsai and Kolsly on the fracture strength 
of glass plates as affected by indentor diameter, by 
observing the anologous nature of both the problems. 
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The present f omul at ion gives an alternative approach 
of accurately characterising the size dependent scatter 
in strength of materials, although leading to mathe- 
matical complexity in computation for materials with 
complicated size-mean strength relations. 

Idealizing cement concrete as a material that 
fits the weakest link concept, an attempt is made to 
characterize the scatter in direct tensile strength hy 
distribution functions using the results obtained in 
the thesis. The experimental results of Kadlecek and 
Spetla have been used in this study. 

As an application of statistical aspects of 
strength in materials testing and specifying strength, 
the minimum number of test specimens required to be 
tested to nredict the mean strength with a given relia- 
bility as affected by specimen size, is investigated 
for cement concrete. Since the coefficient of variation 
of strength decreases with increase in specimen size, 
it is shown that more specimens of smaller size are to 
be tested for a given reliability in predicting mean 
strength, than for large sizes. A study of optimum 
specimen size and sample number so that total cost of 
testing is minimum is also carried out. 

A unified theory is -proposed which correlates the 



size effects on strength, stiffness, fracture mode 
transition and nonlinearity in materials like concrete 
and fiber reinforced mortars in which fracture initiation 
and total failure are not identical (due to the presence 
of microcracking) . Consequently, the weakest link 
concept is not applicable. The concept used is, as the 
applied stress level increases, microcracks are developed 
which cause a progressive breakdown in internal structure 
resulting in the decrease of stiffness. To quantify the 
phenomenon, the probability of failure of the material is 
related to the probable number of microcracks in the 
material. Furthermore, to estimate the loss in stiffness 
i'dth increasing stress, an approximate expression for the 
elastic modulus of a composite material in which one of 
the phases happens to be voids is related to the micro- 
crack volume through a semi-empirical constant. Since 
the probability of failure is related to the applied 
stress level as well as the size, the proposed theory 
yields a size and stress dependent value of the effective 
Young’s modulus, which reflects the nonlinearity in 
stress strain behaviour and the size effect. The proposed 
theory also brings to light a new phenomenon viz, the size 
effect on stiffness and interprets the size effect on 
ductile-brittle transition from a different view point 
from that of Glucklich. The predictions of the theory, 
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i,e. the lower stiffness and lower ductility associated 
with specimens of larger size are shown to he qualita- 
tively in agreement with observations in tests conducted 
on fibre reinforced cement mortar briquettes. 

The fracture resistant design with emphasis on 
material choice, of a pressure vessel in the presence 
of a part-through crack is attempted. The structural 
requirement is that the vessel should withstand the 
internal pressure cycles. The service life is estimated 
by considering the crack growth due to internal pressure 
cycles. In the problem solved, six different materials 
are considered and the material choice is made so as to 
maximize the utility which is defined as the net 
monetary return considering service life, initial cost 
and cost incurred in service because of the weight. 

Problems for further investigations are suggested 
as relevant to those discussed within the scope of the 
thesis. 
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NOTATIONS 


A 

A, 


3/ -j j cl<p ^ « 

a cr 

a final 

a initial 

B 

b 

C 


'T 




= cross sectional area; 

= constant relating fracture stress to crack 
length; 

= cross sectional area of a 6 in. cube; 

= constant in Bolotin’s size - mean strength 
relation; also depth of part - through 
crack; 

a n+1 = roots of the equation V D(Y) = 0; 

— critical depth of part - through crack; 

= final crack depth (critical); 

= initial crack depth (critical); 

= stress dependent risk of rupture; 

= constant in Bolotin's size - mean strength 
relation; 

= crushing stress; also material cost of 
pressure vessel in Rs. per inch length; 

=* total cost of testing in Rs. 

* coefficient of variation; 

* half 'crank length; also semimajor axis 
of part - through crack; 

= crack velocity; 


c^- * half final length of major axis of Crack 

final (critical); 

c initia3 * half initial length of major axis of 
crack (critical); 

c, *= material constants in the distribution 

^ function of Greene; 
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c 

0 

= 

effective fractional void volume; 



constant in size - mean strength relation; 

c 2 

=5 

constant in the extreme value distri- 
bution; 

c 

= 

effective fractional microcrack or void 
volum e ; 

C 1 

=s 

effective fractional microcrack or void 
volume in a specimen of size Y^ ; 

? 2 

sr 

effective fractional microcrack or void 
volume in a soeciraen of size V 0 ; 

D 

=E 

diameter of the pressure vessel; 

D(Y) 


polinomial in Y of degree 'n' ; 

S(7) 

=r 

YD(V); 

B*(.) 

=5 

derivative of the function 5 evaluated 
at (.); 

a 

=S 

lateral dimension of a specimen, diameter; 

d 6 

- 

lateral dimension of a 6 in. cube; 

E 


Young’s modulus; also maximum percent- 
age error of tbe sample average; 

* 

E 

=S 

effective You nr’ s modulus; 

E(.) 

=5 

double exponential distribution of the 
argument (.■) in the theory of Ease; 

sy( r 1 ) 

3S 

cumulative probability of failure up to 
a stress level of a specimen of size Y; 


= 

fracture stress; 

*1 

= 

distribution of smallest values; 

f(x) 

=s 

probability density function of flaw 
length x; 
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) = probability that fracture occurs for an 

applied stress level lying between 
0 and a + cb ; 

J(o ) = function of 0 increasing with 0 and 

satisfying suitable defined conditions; 


■p 

J " o 

f -j > f 2 > • • • 1 ^ 





C- 


G(o' ) 

G c 

G i 

G* 

H( * ) 


h 


= applied force; 

» functions characterizing various 
mechanisms of energy dissipation; 

= functions of 9 characterizing crack 
tip stress field; 

= energy release rate; also shear modulus 
of solid phase; 

= the integral j ) da ; 

= critical strain energy release rate; 

= shear modulus of the i th phase; 

= effective shear modulus; 

= Heaviside unit step function of 
argument ( . ) ; 

= height of specimen; 


h. 


m 


K 

K- 




K 


K. 


= chosen crack length in the theory of 
fisher and Hollomon; 

OO 

- the integral / exp (~Z m ) dZ; 

o 

- material parameter in the size-mean 
strength relationship ; 

= fracture toughtiess; 

= stress intensity factors in the respective 
modes of crack propagation; 

= effective Bulk modulus; 

= cost. of unit volume of material in Rs.; 
also^ return in Rs. from one cycle of 
service of the pressure vessel'; 
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K 


k 


o 


1: 


1 


irk) 


m 


N 


N o, 

N( Y) 
*1 
n 2 

ft 

n 

n(ff ) 


= cost of testing a single specimen in 
Rs, ; also cost incurred due to weight 
in Rs. per cycle of service; 

= length, of pressure vessel; 

= constant in size - mean strength 
relationship ; 

* constant denoting the mean volume of a ' 
single microcrack; 

= scmi-empirical constant relating the 

microcrack volume to effective microcrack 
or void volume; 

= modified stress intensity factor; 

= Inverse Laplace Transform of function (.) 

~ spring constant in compliance method 
of K determination; 

= material parameter in tTeibull’s distri- 
bution function; 

= number oC test specimens; also number 
of phases in a composite material; also 
number of volume elements; 

= service life in cycles; 

= polinomial in V of degree n-1 or less; 

= total number of macroscopic volume 
elements in a specimen of size ; 

= total number of macroscopic volume 
elements in a specimen of size Y 2 5 

= total number of test specimens; 

= number of defects or flaws per unit 
volume; also material parameter; 

= number of elements ruptured during 
loading upto a stress level a ; 



function of stress level ; denoting 
the probability of failure of an element; 

number of specimens ruptured up to a 
stress level a ; 

sample "size; also number of elements 
that get ruptured during loading upto 
a stress level a in a specimen of 
size ; 

number of elements that get ruptured 
during loading upto a stress level 
in a specimen of size VgS 

number of macroscopic elements in the 
material per unit volume; 

prism strength; 

density function denoting the probability 
of the variable (.) talcing a value 
over an infinitismal interval; 

strength of a 6 in. cube; 

probability that a specimen containing 
n flaws has a strength lying between 
(3 and P+dP as a function of a; 

internal pressure, psi; 

crushing load; also flaw shape parameter 
defined in terms of P - and § ; 

critical flaw shape parameter; 

fracture stress; 

constant in Bolotin’s size - mean 
strength relation; 

constant in Bolotin's size - mean 
strength relation; 

distance of point of interest from 
origin in crack tip stress analysis; 



= cumulative probability of failure in 
‘'Veibull's distribution; 

= standard deviation of strength; 

= probability of failure of an element for a 
stress level between 0 and a ; 

= nondimensionalized fracture stress in the 
theory of Fisher and Hoi lam on; 

= nondimensionalized fracture stress; 

= principal stresses along the axes X,Y,Z; 

= strength of material in the presence of a 
flaw; 

= fracture stress in the presence of a crack; 

= fracture stress of a specimen "'ontaining 
n flaws; 

= strength of a material in the absence of 
any flaw; 

= model strength of a specimen in the presence 
of flaw; 

= standard deviation of flaw strength; also 
parameter in Laplace transformation; 

= constant in Bolotin’s equation; 

= mean value of upper yield point stress; 

' = lower bound on fracture stress; 

= surface energy per unit area; 

= time; also denotes function *(u); also 
parameter »t' in the students' * t * distri- 
bution, and also wall thickness of pressure 
vessel; 

= strain energy; 

= net monetary Return in Rs. 



u 


V 

"^1 * *^ 2 * ^ ^ 
TT 

v opt 


V 


V 


o 


V 


v failure 


v 2 


= variable of integration; 

= specimen size volume; 

= volumes of specimens of various sizes; 

= optimum specimen volume; 

= aribitrary volume for nondimensionalizing; 
= volume fraction of voids; 

= volume of microcracks; 

= effective volume of microcracks or voids; 

= microcrack volume at failure; 

= volume fraction of the i th phase; 

= effective volume of microcracks or voids 
in a specimen of size 7^ ; 

= effective volume microcracks or voids in 
a specimen of size V 2 ; 


7 = surface energy; 

w = weight of pressure vessel in lbs, per inch 

length; 

X = random variable; also a coordinate in 

rectangular coordinate system; 

x = flaw strength; 

x Q = lower bound on the value of random variable 

x, = minimum value of the random variable; also 

flaw size; 

# 

x = most probable flaw size; 

X,Y, Z = orthogonal coordinates; 

a = constant in extreme value distribution; 

also denotes the ratio of lateral stress 
to axial stress in a triaxial stress field; 

a n = limit to which f t* tends as u tends to 


oo 
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constant denoting decrease in material 
strength with increase in flaw size; also 
constant defined in tezms of effective 
Poisson's ratio v * ; 

non dimens ionali zed stress along * Z' axis; 
also crack propagation factor; 

constant defined in teims of effective 
Poisson's ratio v* : 

ratio of applied stress to yield stress; 
Gamma function of argument (.); 
space average dilatation; 
i th elemental volume; 

any parameter (like, water - cement ratio) 
on which energy dissipation depends; 

strain at necking instability; 

space average strain tensor; 

nondimensionalized quantity as defined in 
the theory of Ease, denotes the deviation 
of strength, from modal strength; 

angular coordinate of the point of interest 
in crack tip stress analysis; 

constant in the exponential distribution; 

effective lame's parameter; 

mean flaw strength in the theory of 
Prenkel and Kontorova; 

effective lame's parameter; 

Poisson's ratio; 

effective Poisson's ratio; 

dummy variable of integration; 
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% 


a 


a L 

a m 




3>(a ) 


$ (a ) 



= function defined in terms of nondimension- 
alized fracture stress in the theory of 
Fisher and Hollomon; also denotes infinite 
product ; 

= applied stress; 

= nomimal peak stress; 

= stress tensor; 

= lowest mean strength of material; 

= modal strength; 

= ultimate stress; 

- material parameter in 'Weibull's distribution; 
= yield stress; 

= material parameter in Teibull 1 s distribution; 
= mean failure stress; 

= mean failure stress in bending; 

= space average stress tensor; 

= mean failure stress in tension; 

= upper yield point stress in flexure; 

= mean strength of specimen of size ; 

= mean strength of specimen of size 

= Griffith’s fracture stress; 

= probability that fracture does not occur 
for a stress lying between 0 and a ; 

= probability that fracture occurs for a 
stress lying betv/een 0 and cr ; 

= probability density function of flaw 
strength x; 
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- cumulative probability of failure of a 
specimen of size Yy 

- cumulative probability of failure of a 
rpecimen of size 1 0 ; 

- elliptic integral of the second kind; 



function defined in terms of material 
-oarameter a ; 


{<j> (o)d,-}AV = probability that a potential flaw exists 
in an elemental volume aV and is of 
suck size that fracture occurs when 
the applied stress is between a and 
a + dr ; 
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function defined so that f r (u) = 



CHAPTER ORE 


INTRODUCTION 


1 . 1 GENERA1 

The search for improved knowledge of mechanics of 
materials and new design philosophies appears to he a 
continuous process. However, the stringent needs of 
modern technology have emphasized on the need for refined 
understanding of mechanics of materials in order to 
reduce the gap between apparent strength and theoretical 
strength (based on atomic structure and nature of bond), 
new design methodologies, and development of new materials. 
The basic feature that motivated study into all the above 
areas of investigation, is essentially a tendency to move 
closer to reality with minimum of empiricism and uncer- 
tainity. While the earlier attempts of engineering design 
are motivated to?/ards use of already existing materials, 
the current design methods focus on the need for exami- 
nation of material - structure interaction during the 
design process (1,2, 3,4) as well as the design of new 
materials. In addition, the conventional approaches of 
engineering design do not reflect explicit measures of 
design efficiency like optimality with respect to weight 
or cost. However, the stringent demands of economy, of 
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late, require a comparative study of alternative design 
solutions. Consequently, the importance of several 
optimality criteria and cost-effectiveness models has 
heen recognized (5,6, 7,8). 

The search for rationality in the concepts of 
structural safety (9,10,11,12,13) and the observation 
of the randomness in material properties and loads have 
paved the way for considering the reliability as a measure 
of safety. As a consequence, several attempts are being 
mode to examine the probabilistic aspects of strength 
of materials and their testing. Of related interest 
to the aforementioned areas is the phenomenon of size 
effects on material behaviour and associated scatter in 
properties of materials. The importance of size effects 
on material behaviour becomes obvious in any attempt 
to correlate the observations on strength and failure 
mode noticed from simple tests to those in a full scale 
structure. While the earlier studies in scientific 
design processes did not warrant a thorough examination 
of size effects on material behaviour, these cannot 
be neglected in any rational investigation. This is 
especially true of ceramic type materials (for eg: 
clay, cements, gypsum, concretes, porcelain etc.) which 
show brittle fracture characteristics and strength as 
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a function of size. This can be attributed to the 
presence of flairs or internal cracks during their 
manufacturing process. Several studies on the size 
effects on strength of materials have brought to light 
the importance of defects and flaws in materials and 
have paved the way for new disciplines like fracture 
mechanics. As a result of the application of the 
concepts of fracture mechanics several phenomena like 
brittle fracture, fatigue and creep rupture (14,15,16) 
can be reasonably interpreted. In addition, fracture 
resistant structural design can be more rationally 
taclcLed with the use of concepts pf fracture mechanics, 
Svent hough the achievements hitherto, in the aforemen- 
tioned areas like material - structure design, rational 
concepts of safety, statistical aspects of strength of 
materials, standardization of testing methods etc, are 
quite promising, further investigations are essential 
for a better understanding of their nature and application. 
The present thesis is an attempt to investigate some 
of the aforementioned areas and to synthesize the 
current state of knowledge with the current investigation 
on brittle fracture statistics. 

In the following sections of this chapter, the 
salient results of various studies into the areas of 
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fracture mechanics, size effects on strength and 

« 

behaviour of materials and fracture resistant design 
are reviewed while the statistical strength theories 
are critically reviewed in Chapter 2. The presentation 
of the review provides a basis end motivation for 
specific problems that are dealt with in the subsequent 
chapters of the thesis. 

The effect of specimen size on strength and 
mechanical behaviour of materials is of importance in 
engineering design, specification of strength of materi- 
als as well as in correlating observations from tests 
on models to those in the prototype. Some materials 
of engineering interest are more flaw sensitive than 
others, fracture mechanics provides a basis for the 
study of flaw sensitive materials. Event hough flaws 
in materials may be randomly distributed and inter- 
action exists between various flaws, in fracture 
mechanics, only tbe equivalent potential flaw or 
critical flaw and its effect on strength is usually 
studied. Qualitatively, the larger the specimen size, 
greater is the probability of having an effective 
critical flaw in the material and consequently lower 
the mean strength. Consequently the topic of fracture 
mechanics provides the starting point for various 



5 


statistical strength theories. In the following sections 
the basic concepts of fracture mechanics and its appli- 
cation to engineering design are briefly discussed, 

1,2 FRACTURE MECHANICS 

The need for a proper recognition of the presence 
of orach like defects or flaws in materials is felt in 
the wake of several studies on strength of materials 
as well as sudden catastrophic fractures In several 
structures (Refer Sections 1,2,1 and 1,3 of the thesis). 
The mechanism of fracture of materials in the presence 
of crack like defects is termed as fracture mechanics. 

The subject provides a ba,sis for several statistical 
strength theories and fracture resistant design methods, 
as will be discussed later. Vv 

The origin of the evolution of the subject of 
fracture mechanics can be traced bank to the work of 
Griffith (17). Griffith reasoned out that the material 
strength at rupture could be significantly affected by 
the presence of crack like defects in a material. Making 
use of the elastic stress field in an infinite tension 
panel in the presence of an elliptical crack (Fig. 1.1) 
as given by Inglis (18) and using the energy balance 
considerations at the onset of crack propagation, 

Griffith advanced the condition that Unstable crack 
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propagation occurs when the elastic strain energy 
release rate with respect to crack length from the 
medium is at least equal to that of the surface energy 
absorbed by the newly formed crack surface* This can 
be written in the following form, 


cnj > aw 

Be - 0 c 


r, ( 1 . 1 ) 


where IJ * the strain energy 

W = the surface energy 

and 2c = the crack length 

Eqn. 1.1 on substitution of the relevant quantities for 
linear elastic case takes the form 


* 


a 



*r (l t 2) 


where 


and 


a = the failure stress in the presence of arack 
of length 2c 
E * Young’ s modulus 
T « the surface energy per unit area 


Eqn. 1,1 means that fracture occurs when the term on 
the left hand side of Eqn. 1.1 attains a particular 
value called ’Critical Strain Energy Release Rate' 
denoted by G- c , Griffith's energy balance criterion for 
fracture as given by Eqn. 1.1 is applicable only to 
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ideally brittle materials which, behave elastically, m 
which the strain energy is completely reversibly stored 
upto the onset of fracture and no other energy dissi- 
pating mechanism other than surface energy is present* 
It is suggested by Irwin (19) and Orowan (20) that the 
range of applicability of Griffith* s criterion could 
be extended by considering the energy dissipated m 
the form of plastic deformation along with the surface 
energy. As remarked by Narayan Swamy (21), a very 
general form of energy balance criterion for crack 
propagation can be written m the form 


AU — A{f*|(c)+ fg( c , cr , t,S)+ f^(c 


2 2 

t ,E 


,« 2 ) 


- f 4 (c) - f 5 (o 2 , 0 2 )> 


.. (1.3) 


where 



= surface energy, function of crack length c 
- time dependent irrecoverable deformation, 
function of crack length c r stress level 
o»time t and other parameters 6 
= kinetic energy of crack propagation 
dependent on 5 , the crack velocity 
= energy at stress concentrations 
= energy of the applied stress field. 
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Inspite of the fact that Griffith’s work "brought to 
light the significant effect of the presence of for- 
tuitious flaws on material strength, it remained purely 
to he of theoretical value till the real engineering 
theory of fracture mechanics was developed hy Irwin as 
outlined below (22). The lack of motivation to seek 
the applicability of Griffith’s theory to engineering 
problems and design, appears to be due to the fact 
that the quantity, critical strain energy release rate 
G_, lacks a direct physical meaning and also that the 
surface energy cannot easily be measured. 

Before going into work of Irwin and the evolution 
of the engineering fracture mechanics, it is necessary 
to mention the contribution of Sneddon (23). Sneddon 
considered the elastic analysis and the related Griffith 
problem of a penny shaped crack and helped to disclose 
the fact that the elastic stress field m the neighbour- 
hood of a crack can be represented by the stress 
components cr . . , 

cl 

'll - J uq ( 0 ) +k 2 < 0 > +k 3 f i 3 

+ ... other nonsingular terms .. (1.4) 

where = stress intensity factors of the 

dimension psi /“in 
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Q — coordinates characterizing the location 

of the point at which stress field is 
sought (Fig. 1*2) 

1 2 t 

f , f .,f . -= scalar function of 0 . 

1D’"1D ID t 

Enn. 1.-! reflects res oily the so called crack tip 
singularity -ad K 1 , IE,, are consequently known as 
stress intensity factors which are dependent on the 
Young's modulus end crack length. Evaluation of , 

K 2 , m any particular problem is adequate to chara- 
cterize the stress field around the crack tip as far as 
fracture problems are concerned, Sneddon's work not 
only stimulated active research m elastic analysis of 
crack problems by identifying them to fall in a special 
class of boundary value problems, but also paved the way 
for the origin of several physically meaningful results 
m terms of stress intensity factors. An excellent 
account of bhe crack problems in elasticity is given 
by Sneddon and Lowengrub ( 24 ), fhe effect of crack tip 
cohesive forces on the stress intensity factors is 
studied by Barenblatt (26). 

Eoll owing the works of Griffith (17) and 
Sneddon (23), Irwin ( 22 ) brought out the physical 
significance of the stress intensity factors. Irwin 
related these to the concept of critical strain - 



11 


rate, as follows. Considering the case of Griffith’s 
problem, Irwin argued that, by superimposing a self 
equilibrating tensile stress on a free crack surface 
over a length interval (starting from the crack tip), 
the crack can be closed on the same length Interval 
and the force required to do so is the same as the 
strain energy release rate. In other words, the strain 
energy release rate Is the same as ihe crack extension 
force and there exists a critical value of the stress 
intensity factor K called Fracture Toughness corres- 
ponding to & c * A quantitative evaluation of the above 
arguments leads to the relationships 

K q = (Plane stro.cs) 

~Tg c ( 1 . 5 ) 

= — -w (Plane strain) 

V 1-V 

where v is the Poisson's ratio. 

The most general mode of crack propagation in a 
three dimensional medium can be seperated into thre*e 
independent modes of crack propagation (Pig. 1.2), 
namely the tensile mode, the In-plane shear mode 
and the antiplane shear mode. Under general loading 
conditions and crack geometries, the three stress 
intensity factors K.,, Z 2 and are adequate to 
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characterize the :hree nodes of crack propagation. A 
comprehensive catalogue of the values of stress intensity 
factors corresponding to various crack geometries and 
load conditions is compiled hy Paris and Sih (25). The 
fracture toughness of a material can he determined either 
hy employing cracked test specimens (27) or hy the 
compliance method (28,29). In the former, a specimen 
with known crack geometry is loaded and the load and 
crock length at failure are noted and the fracture 
toughness is obtained through the corresponding stress 
intensity factor. In the compliance method the spring 
constant M of a cracked specimen is noted corresponding 
to various values of crack length c ^nd the strain 
energy release rate is given by 


Ct 


( 1 , 6 ) 


where G- = strain energy release rate 

U = strain energy 

f Q = applied force 

knowing G- from Eqn. 1.6 at the onset of fracture, Gk 

7 c 

is obtained and K from Eqn. 1 . 5 . 

V-/ ' 
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1.2.1 Role of Fractur e Mechani c s in Engineering Design 

Several studies on structural elements like 
pressure vessels, ship plates, turbine rotors etc,, 
indicated that tvpical brittle failures occuring at 
operational or service loads are originated by small 
cracks or crack like flaws (63,64,65). A check for 
fracture resistance m structural design assumes gfeater 
importance while using high strength materials since 
fracture toughness usually decreases with increasing 
strength and as such the design conditions are more 
rationally specified m terms of fracture stress rather 
than yield or ultimate stress. Of the three approaches 
to fracture resistant design viz, the transition temper- 
ature approach (67,68), the application of notch 
stresses and the fracture mechanics approach, the last 
mentioned one is more general in its application. 

The fracture mechanics approach to design of 
structures is based on the criteria and conditions under 
which self supported crack propagation can take place 
as dependent on the material properties, presence of 
crack like imperfections, size and shape of the body 
and other environmental conditions. G-iven a design 
problem and the possible materials, use of fracture 
mechanics approach m principle consists of the 



following steps. First the fracture toughness of the 
various materials at various temperatures is determined, 
the temperature ranges being relevant to the design 
problem. The central problem of the solution is in 
estimating the most likely critical crack sizes, shapes 
and orientation, that can occur during the service life 
of the structure as affected by the sustained and fati- 
gue loads. While the crack size estimates can be made 
depending ~>n the accuracy of the nondestructive testing 
equipment available, fairly reliable information of 
engineering value can be obtained from full scale 
testing and a study of the earlier case studies. 

The next phase of the design process is an analy- 
sis of the structure for a crack size-strength relation- 
ship at various temperatures. Since the stress intensity' 
factors as relevant to various crack shapes and loading 
conditions of practical interest are available ( 25 ) 
and various numerical techniques are possible for their 
ready evaluation, the aforesaid analysis poses no 
problem. With this information as well as economic 
and other design values, the material that suits best 
is chosen keeping an estimate of the most probable flaw 
size. The aforementioned procedure is essentially 
cursory in exposition and various design problems call 
for special techniques in adopting the basic principle 



outlined. A detailed account of the various design 
practices as applied to different structures can be 
found in the work edited by Leibowitz (64). An inte- 
grated fracture design approach with a measure of 
optimality is discussed in Chapter 9 of the thesis. 

1,3 SIZE EFFECTS CBT STRENGTH OF MATERIALS 

Classical concepts of flow and fracture of materi- 
als based on various postulates like limiting state of 
stress or limiting state of strain criteria are based 
on certain idealizations of material behaviour and as 
such are inadequate in rationally visualizing various 
anomolies normally experienced with failure behaviour 
of materials. This is because the flow and fracture 
of materials in its origin is a highly localized pheno r 
menon basically initiated and propagated by various 
Inherent defects or flaws or structural irregularities 
present in almost all materials. Defects could be both 
on macroscopic and microscopic scale. Lattice imper- 
fections in crystals, dislocations and grain boundaries 
in metals, surface flaws in glass, microcracks and 
inclusions of foreign materials (like that of aggregates 
in concrete and minerals in rocks) are typical examples 
of these strength impairing defects. Eventhough the 
exact origins of imperfections are not clearly known 
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yet, these could he considered as inherent material 
characteristics which are originated during the process 
of material manufacture and subsequent treatment given. 

By their nature, imperfections are random in number, 
orientation and distribution throughout the material. 

The classical approaches to flow and fracture- attempt 
to examine the failure behaviour only on the average 
and as such are mathematically simple. However, such 
attempts cannot explain the phenomena like size effects 
on strength, scatter associated with testing of 
materials for static and fatigue conditions etc. 
Furthermore, the scatter associated with testing of 
materials, specimen size effects etc. are of significant 
engineering importance. This can be visualized easily 
by noticing that brittle materials display strength 
variations ranging from about 30% to 30C% of their 
mean strength (14-), even under apparently identical 
test conditions. Consequently, characterizing such 
materials by a single value called mean strength may not 
always be justifyable in terms of engineering application. 
Recognizing that materials are flaw sensitive and that 
flaw occurance is not a deterministic phenomenon, one 
is prompted to the application of mathematical theories 
of probability and statistics. Before going further 
into this aspect of the problem in the thesis, an attempt 
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is made in what follows to review and summarise several 
investigations of size effects on strength in different 
materials. The statistical aspects of the various 
studies are critically reviewed in Chapter 2 , 

1 . 3*1 Cast Iron and Mild Steel 

McPherran (30) experimentally investigated the 
effect of increased cross section size on the tensile 
strength of cast iron. Three types of the material were 
employed, in each case the diameter of the test speci- 
men was the variable. The diameters were varied from 
1.25 in to 4 in and intotal specimens of six sizes 
were tested. It was noticed that increased specimen 
size caused a decrease in the mean tensile strength. 
Campbell( 3 l) conducted tension tests on cast iron rods 
of 15 in. length and of diameters 0.25 in., 0.375 in., 
0.625 in., 0.875 in., and 1.2 in. The specimens were 
prepared from nine different heats and the results 
indicated a decrease in tensile strength with increased 
specimen size. Schneidewind and Hoenicke (32), in a 
very comprehensive experimental investigation have 
studied size effect on strength of gray cast iron. 

For assessing the tensile strength of the material 
the authors have employed test hars of sizes 0.5 in., 
0.875 in,, 1.2 in. and 2 in. in diameter and of length 
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18 In. to 24 in. The observed tensile strengths of the 
bars respectively were 52,267* 34,142^ 32,050 and 
22,523 psi. Strength of the some material in compression 
was studied by testing cylinders in compression with 
length to diameter ratio of two, and the strengths were 
respectively 112,342, 112,017, 119,583 and 99,783 psi 
which show lesser variation than that for the case of 
tension. The modulus of rupture of the material *was 
studied by testing the rods of various aforesaid diameters 
on spans of 8 in., 12 in., 18 in. and 24 in. respectively, 
with centre point loading. The corresponding values 
of modulus of rupture were 98 , 950 ; 75,933; 70,783 and 
62,700 psi. All the above investigators have concluded 
that the size effect on the observed behaviour is 
essentially because of the effect of bulk of material 
on the rate of cooling from molten state to room temper- 
ature In the preparation of the test specimens and the 
consequent change in the Internal metallurgical 
structure of the material. 

Davidenkov et.al (33), in their study on size 
effect on brittle strength of steel in liquid air 
atmosphere noticed that the mean strength and dis- 
persion decreased with increase in specimen size. Thrther 
the variation of strength and dispersion was adequately 
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explained "by Feibull’s theory. Richards (34) examined 
the size effect in the tension testing of mild steel 
with emphasis on the upper yield point. In an attempt 
to interpret the existence of upper and lower yield 
points, Richards proposed a model for mild steel, 
which is explained below. He considers that the material 
consists of two components one brittle and the other 
soft or ductile. Under increasing loads, the system 
behaves elastically up to upper yield point at which the 
brittle component suddenly fails e.nd the load is taken 
by the soft component only. As a consequence, Richards 
considers that the phenomenon is governed by the same 
laws as that of brittle fracture, like the occurance 
of random inhomogenties, microcracks etc. The series 
of tests conducted on rods of 1/8 in., 1/2 in. and 1 1/2 
in. diameter with corresponding volume ratio 1:64:1000 
have shown values of upper yield point being 60,3 00, 
55>600 and 53,550 psi. It was found that the relation 


60. 14 

u “l/5870 


.. (1.7) 


where s^ = mean value of the upper yield point stress 
in psi 

V = specimen volume in cubic inches 


fits the experimental data well 
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Richards (35) in a subsequent investigation of size 
effects on yielding of mild steel in flexure corrobo- 
rated his earlier findings. He tested beams of mild 
steel of five different sizes, ten in each size, ill 
the specimens of different sizes were made dimensionally 
similar, the ratio between the corresponding dimensions 
of the smallest and largest specimens being 6.3 while 
that between adjacent being 1.59. The results indicated 
that size effect was present on the upper yield point 
of mild steel and the following relationship 

? ‘ - -22222 .. ( 1>8) 

t 11.7 

where = the upper yield point in flexure in psi 

3 

Y = the specimen volume in . 

adequately describes the variation of mean value of 
upper yield point with specimen volume. 

1.3.2 Concrete and Rock 

One of the earliest investigators of the size 
effect on compressive strength of concrete was by 
Cornerman (36). Based on tests conducted on cylinders 
of diameters ranging from 1.5 in to 10 in. with length 
to diameter ratio 2, he noticed a decrease in strength 
with increase in size. Reagel and Willis (37) in their 
study of size effect on modulus of rupture of concrete 



21 


beams noticed that the variation in the modulus of 
rupture was negligible with increase in length or width 
of the beams. Modulus of rupture was found to decrease 
with increa.se in depth; the decrease in modulus of 
rupture was about 1.9% for each increase in depth by an in. 
Kellerma" (38) in his investigation noticed that, the 
flexural strength of concrete decreased with increase 
in span as well as the cross sectional area. Tucker (39), 
in his analysis of the size effects on modulus of 
rupture, suggests that the weakest link concept adequately 
interprets the decrease in modulus of rupture with 
increase in span length and depth while a parallel element 
model is essential to interpret the fact that increase 
in beam width has neglible effect on modulus of rupture. 
Johnson (40), in his study of scaled down concrete f inds thatfrr 
a given scale of mix the mean compressive and tensile 
strength decreases with increase in specimen size, the 
variability being more at smaller sizes. Raftendran (41) 
tested concrete cubes of sizes 4 in. , 6 in., 8 in. and 
10 in. in compression. In each case twelve specimens 
were tested and it is noticed that the strength decreased 
with increased specimen size. The test results indicated 
that taking the strength of a 4 in. cube as 100% the 
percentage strength of 6 in. cube was 95.62%, that ef 
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8 in. cube 92.14% , that of 10 in. cube was 87.84% . 

The trend indicates that the strength decreases by 
about 4% for each 2 in. increase in the size starting 
from 4 in. cube. In an attempt to find an overall 
relationship between the strength of concrete and its 
size and shape, Neville (42) identified three parameters 
viz, d the maximum lateral dimension, V the volume and 
h/d the height to lateral dimension ratio. Making 
use of the test results of several investigators on 
cubes, cylinders and prisms, Neville proposed the 
foil owing rel ati onship s : 

— = 0.56 + 0.697 .. (1.9) 

P 6 (& - *) 


, A , 0.4525 
0.88 78 (f-) 

A 6 


( 1 . 10 ) 


where P = the strength of the specimen psi 

Pg = the strength of a 6 in cube psi 

2 

A * the cross sectional area in 

2 

A * areaof a 6 in. cube in 

6 

v v 

Kadlecek and Spetla (43) have made an extensive study 
of the size and shape effects on direct tensile strength 
of concrete. They have tested cylinders and prisms of 
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different volumes and of constant length, to lateral 
dimension ratios. The tests indicated a decrease in 
tensile strength and the associated scatter with increase 
in specimen size. The analytical results based on the 
data of this experimental study are discussed more in 
detail subsequently in the thesis. Jaeger and Cook (43) 
discuss available data on rocks and conclude that the 
tensile strength of rocks is a highly variable property,-' 
and is influenced by specimen size more than any other 
property of rocks, 

Skinner (44) studied the effect of specimen size 
on compressive strength of Anhydrite and noticed that 
the moan strength and standard deviation decrease, with 
increase in test specimen size. 


1.3.3 Class. Alumina and Porcelain - Ceramic Materials 
Griffith (17) is the earliest investigator to 
have studied the size effect on strength of glass 
fibres in a formal manner. He reported that, within 
the limits of experimental error, the equation 


= 22,400 £&±-§ 


22,400 -^8600 


.. ( 1 . 11 ) 

adequately reflects the effect of diameter d (inches) 
on failure stress (psi) H 0 of the material. Eqn. 1,11 
is valid within the range of diameters that have been 
considered in the tests viz, 0.001 in. to 4.20 in. 

Size effect on tensile strength of plexiglass (although 
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not a ceramic material) , Teas teen studied "by Burelli 
and Paries (75) using dogbone specimens in tending and 
tension and the applicatility of the relation 

F q = c Y -lc . . (1.12) 

where F = the failure stress 
o 

c and 1c = constants 

7 = the specimen volume 

Greene (46) studied size effect on strength of glass t y 
testing circular rods in flexure. The investigation 
showed the size effect as determined ty the nature of 
surface polishing and the surface area of test specimens. 
The results of Durelli and Parks as well as those of 
Greene will further he discussed in the light of 
statistical strength theories. 

Weil, Bortz and Firestone (47) studied the effects 
of prior thermal history, specimen finish, test temper- 
atures, environment and specimen size on the strength 
of Alumina, lest specimens of doghone shape of 
volumes 0.0117 cu.in. , 0.0469 cu in. and 0.0977 cu in. 
have been used in the test programme. It is noticed 
that the average tensile strengths corresponding to the 

three sizes are 24.4 x 10“ 3 psi, 27x10~ 5 psi and 

—3 

19.1 x 10 psi; while the corresponding standard 
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deviations are 5.7 x 10 5.6 xIO ^ and 2.7 x 10 ^ psi. 

Size effect on bending strength of porcelain has 
been studied by Weibull (48) and Salmassy et.al (49,50). 
Essentially the material displays decrease in strength 
with increa.se in size. Since these studies are of 
significant importance in the development of new 
statistical strength theories, these studies are 
discussed subsequently, in detail. 

1.5.4 Other Materia ls 

The investigations of Comben (51) on timber beams 
with varying depths from 8 cm to 0.5 cm indicated size 
effect on modulus of rupture. The variation in the value 
of modulus of rupture was between about 10000 psi to 
IJOOO psi corresponding to the range of sizes of beams 
tested. 

Baratta and Driscoll (52), during an attempt to 
develop a new tension testing device that minimizes the 
inherent M parasitic” bending stresses, have investigated 
the size effects on strength of graphite. Specimens 
of different volumes ranging from 0.0547 cu in. to 
0.4552 cu in. were tested and the corresponding average 
tensile strengths varied from 9600 psi to 8700 psi. 
llso, the coefficients of variation were between 7 to 10^- . 
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Size effects on crushing strength of coal has teen 
investigated hy Evans and Pomeroy (53). The authors 
have tested coal cutes of sides 0.25 in., 0.5 in., 1 in. 
and 2 in,, the number of cutes tested teing 262, 164, 62 
and 25 respectively. The mean crushing strengths 
correspondingly were 3700 psi, 2800 psi, 2380 psi and 
1820 psi. The standard deviations were 1530 psi, 840 psi 
and 710 psi and 750 psi respectively. Based on the 
atove observations Evans and Pomeroy recommonded the 
applicability of the following formulae 

c a d -0.32± 0.02 .. (1.13) 

q a d 1.68± 0.02 ‘ .. (1.14) 

where C = the crushing stress psi 

d =* the side of the cute, inches 
and q = the crushing load Its. 

Studies on gypsum mortar for modulus of rutpfure 
and compressive strengths on cylinders, by White and 
Satnis (56) have indicated the size effect on strength. 
The decrease in strength with increase in size is more 
pronounced in the case of specimens that were cured in 
an unsealed condition rather than those in sealed 
condition^ This feature was partly attributed to the 
variation in the material as affected by the nature of 
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curing. 

Test results on size effects on strength of stearic 
acid, plaster of paris, cotton threads as reported by 
Oeibull (48) indicate the same trend of decrease in 
strength following increase in specimen size. Some of 
these results are subsequently discussed in the thesis.. 

Shaffer (54) in his study on the growth and 
properties of whiskers discussed the size effects on 
tensile strength of lithium flcuride crystal and 
sapphire whiskers and notices that there is a significant 
decrease in strength with increase in crystal thickness, 
area as well as circumference of the whiskers. 

Some empirical aspects of the size effects in 
various material systems with respect to strength have 
been reviewed in detail by Sabnis and Aroni (55). It 
can be seen from the aforementioned studies that the 
general feature of size effect on strength, is that the 
strength as well as the associated scatter decreases with 
increased specimen size. This aspect of the scatter is 
further examined analytically in detail subsequently. 

1.4 SIZE EFEEC TS OR MEC HANI CAL BEHAVI OUR 

Of wide interest in the mechanics of materials 
is the size effect on strength and mechanical behaviour. 
While the effects of temperature, rate of loading and 
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state of stress on mechanical behaviour have been studied 
to some extent, the size effect is relatively less 
understood. However the size effect is important as 
it is essential in correlating the behaviour of test 
models and prototypes in engineering design in terms of 
performance as well as in establishing safety (or 
reliability) limits. If the occurance of imperfections 
in materials affects strength it is a natural consequence 
to note that the same imperfections might affect other 
aspects of material behaviour such as stiffness, 
ductility etc. 

Schneidewind and Hoenicke (32), in their experi- 
mental investigation on gray cast iron in which they 
tested bars of diameters 0.5 in., 0.875 in., 1.2 in. and 
2.0 in. have noticed that the corresponding Young’s 
moduli E of the material are 19530000, 17570000, 13600000 
and 11980000 psi respectively. The trend indicates that 
the Young’s modulus E decreases with increased specimen 
size. The variation in the value of Young’s modulus 
with specimen size is once again attributed to the 
effect of rate of cooling of the material and conse- 
quent variation in internal structure of the material* 

In materials like concrete it is known that the 
specimen size affects the surface drying and consequently 
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the creep characteristics of the material. Troxell, 
Raphael and Davis (57) in their tests for creep on 
concrete using cylinders of diameters 6, 8 and 10 in. 
with length to diameter ratio 2, have noticed, creep 
ratio decreased with increased specimen size i.e. 
smaller specimens creep more than larger ones. 

Several investigations on fatigue of metals 
indicate that fatigue strength decreases appreciably 
with increase in specimen size. The size effect on 

i 

fatigue is attributed to the occurance of flaws and 
methods of manufacture (15, 58). 

nf great practical importance is the fracture 
mode in materials as affected by specimen size. The 
bulk of the specimen influences the initiation and 
propagation of flaws in the materials and it is known 
that the fracture toughness of materials decreases 
significantly with specimen thickness (59). This is attri- 
buted to the change in state of stress from plane 
stress to plane strain and therefore the extent of 
plastic zone around the crack tip. Qualitatively, the 
zone of deformation will be smaller in the case of plane 
strain and consequently lower fracture toughness (60). 
Decrease in fracture toughness with increase in thick- 
ness is of importance in engineering design as it 
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influences the fracture mode from ductile to "brittle 
as the thickness increases. 

The size effect on nonlinearity in stress strain 
behaviour and ductile brittle transition is relatively 
studied less. Experimental and theoretical studies 
due to G-luclcLich (58,61,62) a.re the few available. A 
formal disctission of these is carried out in detail 
in Chapter 7 in the thesis wherein a theoretical model 
to explain the phenomena is developed. 

1.5 OBJECT AID SCOPE OF THE THESIS 

Any problem in structural design can be considered, 
broadly, to consist of the three choices, i) choice of 
configuration ii) choice of material and iii) choice 
of failure mode. Of basic importance in a design process 
involving the three choices is the knowledge of behaviour 
and properties of materials. Host engineering materials 
are flaw sensitive and the flaw occurance being usua- 
lly c. random phenomenon is size dependent in addition 
to other parameters, such as temperature, state of stress, 
shape of specimen, environment etc. Since the basic 
properties of any material like strength, stiffness 
and mode of fracture are determined from simple tests 
in laboratory on small specimens and are used to 
predict the properties in a prototype, size effects 
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play a significant role in correlating the observations 
from simple tests to those in structures. In materials 
where the scatter in strength is significant because 
of random flaw occurance, (for example: cement concrete) 
specification of strength by a single quantity has no 
meaning and the strength has to be specified by a 
probability distribution function. A knowledge of the 
true strength distribution function is essential in 
rational safety specification by applying probabilistic 
structural design methods. Svent hough, some methods 
like Ang f s (12,13) extended reliability method are so 
formulated that they are insensitive to the exact form of 
distribution functions involved by a synthesis of 
subjective and objective probabilistic and deterministic 
concepts, a knowledge of the true strength distribution 
functions is not unwanted as they enable a more refined 
application of the design method. A consequence of 
scatter in strength of materials, as being affected by 
size of material specimen, is a need to plan a testing 
programme with a decision on the number and size of 
specimens to be adopted in studying the properties of 
a material depending on the state of stress, especially 
for truely brittle materials, 

A knowledge of the mode of fracture of a material, 
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ductile or "brittle is necessary in predicting the 
"behaviour of a prototype structure using the behaviour 
of a model. The phenomenon of size effect on ductile- 
brittle transition was first recognized by Glucklich 
(56,61,62). The importance of this phenomenon in 
predictive testing can hardly be over-emphasized. 
Theoretical and experimental investigations into this 
phenomenon would be helpful in assessing the material 
behaviour as well as in design applications. 

Brittle fracture being catastrophe and sudden 
in nature with a consequent damage, design methods with 
emphasis on material choice (so that a predefined 
utility of the design decision is maximized) need to be 
studied. 

As could be seen from the above discussion, the 
various areas of study outlined would be of theoretical 
as well as of practical importance. In this thesis 
an attempt is made to investigate some problems in the 
areas of statistical aspects of size effects on 
strength and fracture behaviour of concrete-like 
composite materials and fracture resistant design. 

The problems studied in the thesis and the arrangement 
of chapters in the thesis is as follows. 
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In Chapter 2, a critical review of the various 
existing statistical theories of strength of materials 
and the limitations are discussed in detail. It Is 
emphasized that the practice of assuming a priori, the 
form of strength distribution function and adjusting the 
materials parameters to a particular material is 
inadequate in some cases. The need for an alternative 
approach, vis, a method of constructing a strength 
distribution function corresponding to a size-mean 
strength relationship of any material as obtained from 
tests, is indicated. 

In Chapter 3, the formulation of the problem of 
size effects on material strength and scatter is carried 
out using the concept of weakest link theory and the 
random flav T occurance in materials. It is shown that 
the determination of the size dependent strength 
distribution function in materials amounts to the solution 
of a nonlinear integral equation that can be specialized 
depending on the nature of material under consideration, 
i.e. the size-mean strength relation. The motivation 
for the study and the formulation is the work of Tsai 
and Ilolsky (87) on the effect of indentor diameter on 
the impact strength of glass. The correspondance bet- 
ween the problem of Tsai and Kolsky and that of 
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size effects on strength is utilized in formulation 
of the problem. 

In Chapter 4, the method of solution of the 
integral equation on the lines of that of Tsai and 
Kolsky (87) for strength distribution function is 
developed in detail. A closed form solution of the 
equation yielding a, distribution function that can be 
specialized to ^eibull’s (48) or Bolotin's (72) form 
is obtained. A discussion of the possible numerical 
approaches to solve the problem corresponding to more 
general size-mean strength relation is then made. 

Chapter 5, deals with the characterization of 
direct tensile strength of concrete by a statistical 
distribution function. Direct tensile strength of 
concrete has not been examined by earlier investigators 
in the light of statistical strength theories. In this 
Chapter, making use of the experimental results of 

S/ o' 

KadLecek and Spetla (43) and the theoretical results 
of the preceding Chapter 4 of-tlie thesis, material 
parameters in the distribution functions are found by 
trial and error) and a comparison of the theoretical 
and experimental distribution functions is shown ■ 
graphically. 
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In Chapter 6, the effect of specimen size on the 
strength and scatter in materials testing and specifi- 
cation of strength is considered. Making use of 
Students t distribution, the minimum number of test 
specimens required to be tested as affected by size, 
to specify the mean strength with a given maximum error 
and reliability in prediction is discussed. Test 
results of Kadlecek and Spetla (45), are used to illu- 
strate the effect of specimen size on the number of 
test specimens (in direct tension testing of concrete) 
required to be tested for a given error and reliability 
in prediction. It is observed that the number of test 
specimens increases with decrease in specimen size. 

Since the total cost of any testing programme is 
dependent on the quantity of material used as well as 
the number of tests, a study is made of the optimum 
sample and specimen size so that the total cost of 
testing programme is a minimum. 

In Chapter 7, a phenomenological theory is 
developed to interpret size effects on stress-strain 
and fracture behaviour of concrete-like composite materi- 
als in which subcritical crack growth before rapture 
is significant because of the presence of energy 
dissipating mechanisms like microcracking. The theory 
is developed on the concept that at increasing stress 
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levels, the probability of failure of materials 
increases and the probability of failure is related 
to the internal mic roc racking of the materials through 
a semi empirical constant. Using a result for the 
Young’s modulus of a two phase materials, in which 
one of the phases happens to be voids, an expression 
for size and stress dependent Young’s modulus is 
derived. 

In Chapter 8, an evaluation of the predictions 
and implications of the theory developed in Chapter 7 
are discussed. Since the theory correlates probability 
of failure, stress level and Young’s modulus, it is shown 
that size effects and nonlinearity in stress strain 
behaviour are built in the development. The theory is 
shown to predict lower stiffness associated with speci- 
mens of larger size. Considering the energy under the 
stress strain curve as a measure of ductility, it is 
shown that the theory predicts larger specimens display 
smaller ductility. The implications of this theory are 
compared with those of G-lucklich’s theory of strain 
energy and size effect on ductile-brittle transition 
in materials. The predictions and implications of the 
proposed theory viz, lower stiffness and lower ductility 
of specimens of larger size are observed to be quali- 



tatively in agreement with the observations in tests on 
fibre reinforced cement mortar briquettes tested in 
tension. 

In Chapter 9, fracture resistant design of a 
pressure vessel in the presence of a part-through crack 
with emphasis on optimal material choice is attempted. 
After emphasising on the various embrittling factors 
in pressure vessels and the need for design for fracture 
resistance, a specific design problem is considered in 
which the internal pressure and nominal working stress 
level is specified*. The vessel is considered to have a 
part-through crack in the Avail of the vessel. It is 
required to choose a material out of six available 
materials so that the utility of the choice is maximum, 
Thile considering the utility, the material cost, weight 
and service life are taken into consideration. Two 
aspects of the design problem are considered viz, the 
safe life design and fail safe design* 

In Chapter 10, a summary of the work presented 
in the thesis is carried out. ' Further applications 
of the significant results of the thesis are pointed. 
Also potential areas of further study are indicated. 
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CRITICAL REVIEW OE STATISTICAL THEORIES OP 

— s tre^it^ MATERnrr” 


2*1 IrTTRODU CTI OF 

As could "be observed from the investigations of 
various research workers on size effects on strength of 
materials (reviewed in Section 1 . 3 of the thesis) it is 
a common observation that specimens of smaller size have 
higher strength. Further, a general feature of the 
various studies reported, is that scatter in strength is 
characteristic of wide class of materials and that the 
mean strength approaches asymptotically a certain minimum 
value as the size increases. Also the scatter in 
strength values for smaller sizes is larger than that 
for bigger sizes. These features are shown schematically 
in Pig. 2.1. 

The phenomena of size effect and scatter prompted 
several investigators to advance various statistical 
theories of strength of materials. The basic problem 
of statistical theories of strength is to evolve a 
statistical distribution function that adequately 
characterizes the scatter in the strength of materials. 

The problem in effect is to identify admissible forms 
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of distribution functions with suitable parameters 
that reflect accurately the true material behaviour for 
which the distribution function is applied. Also, any 
proposed theory should be on the basis of a realistic 
mechanism of failure. Several theories that are developed 
on these guidelines are reviewed in the following. The 
limitations and merits of various theories are discussed 
subsequently at the end of this chapter. 

2.1.1 Weakest link and Classical Bundle Concepts 

Basically, two distinct approaches are available 
to study the statistical aspects of strength of materials 
and size effect. The approaches are based on the weakest 
link and classical bundle concepts (14). In the weakest 
link concept, the presence of a severe defect anywhere 
is adequate to cause total failure; that is failure of 
total material is same as the failure of any of the 
constituent elements. Consequently in this approach, 
strength In the large is determined by the strength of 
the weakest element present. On the other hand, in the 
classical bundle concept the strength is not determined 
by that of the weakest element alone but is dependent 
also on the strength of the elements in the neighbourhood. 
In this model the material is supposed to be replaced 
by a bundle of parallel fibers and in such a situation 
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the gross strength at failure is influenced by the 
strength of all constituents fibres. These two approaches 
are, in a way, idealizations to make the problem tractable 
and in reality the actual characteristics of materials 
fall in between these two idealizations. Theories based 
on weakest link concept are reviewed in this chapter, 
while those based on bundle concept can be found in the 
works of Daniels ( 69 ) and Druedenthal (14). The weakest 
link concept has been assumed in the development made 
in Chapters 3,4 and 5 , while a theory of composite 
materials which lies in between the weakest link and 
classical bundle concept is used in Chapters 7 and 8 for 
a unified theory of concrete -like materials for fracture 
as related to size of the specimen. 

2.2 WEI BULL 1 S THEORY ADD APPLICATIONS 

2*2,1 Basic Assumptions and Development of the Theory 

Wei bull is one of the earliest proponents of a syste- 
matic study of the statistical aspects of strength of 
materials (48,70, '71). Weibull pointed out the inadequacy 
of characterizing material strength by a single quantity 
as is usually done in deterministic approaches and built 
up a theory in which additional parameters are invoked 
to characterize the strength of a material. Weibull 
employed the weakest link concept in that he considered 
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the failure of the total material is the same as the 
failure of any one of the constituent elements or the 
weakest element. If the material is considered to 
consist of a series of elements of unit volume’ and if 
the probability of failure of the element for a stress 
lying between 0 and o is S , then the probability of 
survival of the element is given by (1-S Q ). If S denotes 
the cumulative probability of failure of a specimen of 
total volume V, then the probability of survival of the 
total specimen is given by 

(i - S) - (1 - S o ) Y .. (2.1) 

or log (1 - S) = V log (1 - S Q ) .. (2.2) 

Weibull defines a function B termed as 'risk of rupture' 
given by 

B = -log (1 - S) = -V log (1 - S 0 ) .. (2.3) 

The risk of rupture of the element is given by 

dB = -log (1 - S 0 )dY .. (2.4) 

It is reasonable to assume that S Q is some function 
n Q (cr ) of stress level a so that 

dB = n Q (a) dV .. (2.5) 


or 


B * / n 0 d) dV 


• ♦ 


( 2 . 6 ) 
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- f n(o ) dV 

and (1 - S) = e 

- f n(o) dV 

so that S = 1 - e 


.. (2,7) 
( 2 . 8 ) 


It can he seen that from the above expression for the 
probability of failure S, knowing n 0 (cr) one can readily 
evaluate S. 

for the case of a simple uniaxial stress field, 
Weibull intuitively suggested the use of the form of 
function n Q (a ) given by 

n o 0O = (f ) m .. (2,9) 

o 

where m and o Q are two additional parameters to chara- 
cterise the strength distribution of the material; m 
is called the flaw density parameter while a Q is 
called the scale effect parameter. 

With this assumed form of n (<*), Eqn. 2,9, S is 
given by 

S * 1 - exp { -V (— ) m } .. (2.10) 

°o 

It can be seen from the above equation that the speci- 
men of volume V has a probability of failure unless 
the stress level is zero. This implies that the material 
can have the lowest strength as zero. However many 
brittle materials do sustain some minimum stress 
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However, Weibull ( 48 ), in his paper expresses a in terns 
of an integral 

I = f e ^ dz Equivalent to Eqn. 2.13 
m ^ 

2 

The variance of strength is given by 

S d ^ = ° Q 2 y“ 2 / m {r (1+2/m) •- r 2 (l+ ~ ) } .. (2.14) 

where the complete gamma function is defined as before- 

_ 2 

The above expressions for cr and indicate 
that as specimen size increases the mean strength and 
variance decrease in agreement with the empirical 
observations of size effect on some materials. 

yeibull further generalized his theory to multidi- 
mensional state of stress and varying stress fields. In 
these extensions, he neglected multiaxial stress inter- 
action (unlike that done- in strength theories for complex 
state of stress) and assumed that only principal tension 
initiates failure. In the present thesis the case of 
uniaxial stress field is considered and as such no 
attempt is made to study the extensions of Weibull for 
other cases. 
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2.2.2 Bolotin. 1 s derivation of Weibullfa Distribution 
Bunc'tion ~*'~ 

Bolotin (72) rederived Wei bull 1 s distribution 
function in a more formal manner by invoking a result 
of the distribution of minimum values. Using the basic 
concepts of the extreme value statistics, Bolotin shows 
that the distribution F^(x) of minimum values x^ of a 
random variable X, for very large sample size n^ , can be 
given by the following 


cc 

E-jCx-j) = 1 - exp {- c 2 n 1 (x 1 -x o ) } for x 1 > x Q 


- 0 


for x. < x 
l-o 


..(2.15) 


where c 2 , a = positive constants 

and x Q = the lower bound on variable X 


Identifying the problem of strength theory of 
materials based on the weakest link concept as the 
same as the study of the minimum strength of a defect 
in a material in the presence of a population of 
defects, Bolotin uses the analogy to write the strength 
distribution function as follows 


) = 1 - exp 

* 0 


{-c 2 nY (R.,-s: o ) a } 


for R^ 
for R.| 


> 


s 


0 


< s o 


..(2.16) 



where 


= cumulative probability ox failure at a 
stress level of a specimen of size V. 
Y = volume of the specimen 
n - number of defects per unit volume 
= fracture stress 

s Q = lower bound on fracture stress. 


Eqn. 2.16 in a nondimensional form is written as 


VV = 1 - exp (^-~) a > •• (2 * 17) 

o c 


where c 2 n is replaced by 


c 2 n = — ~ . . ( 2 . 18 ) 

T . a 
Ys n 
o c 

where. s c = a constant 

The size~mean strength relation corresponding to the 
above distribution function is given by 

y 1/a 

£ * s o + s c ^ r(l+ 1/a) .. (2.19) 

where H = mean strength 


Bolotin remarks that the size-mean strength relation 
for a wide class of materials can be given by a for m 
such as 
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where R . a, b and a are constants. 
o 7 7 

By comparing the two Eqns. 2.19, 2.20 


r ( i+i/a) 


Bolotin shows that 

.. ( 2 . 21 ) 

.. ( 2 . 22 ) 


The coefficient of variation Cy of the distribution 
function Fy(R^) is given by 

V 1/a 

■T 


C Y *- 


b (v £ ) ' * 0 ( a ) 


a + b(y) 


T o. V 0 


.. ( 2 . 23 ) 


where 


4 > 0 ( a ) 


_ i 

\l r 2 (i+ 1/a) 


.. ( 2 . 24 ) 


Gy - coefficient of variation 


It can he noted that the results of Bolotin and Weibull 
are identical except for the method of approach, in that 
Weibull assumed the form of distribution function, 
while Bolotin invoked the result from extreme value 
statistics. 
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2 . 2 , 3 Reduct ion of D ata Suitab le, t o Weibul^s Distribution 
Function and Appli cations 

While there exists no formal basis for the assumed 
form of the distribution function of ,Tr eibull , the remark- 
able feature of the theory is its simplicity far 
application. As has been mentioned already, the basic 
problem of statistical strength theories is to fit 
experimental data into a suitable distribution function. 

If the form of distribution function is fixed a priori 
as is done by Weibull, the material parameters a Q , a u 
and m should be adjusted to have the best statistical 
fit for scatter observed in experimental determination 
of strength of particular material. The cumulative 
probability of failure, according to Eqn. 2.11, is 
given by 

a - a m 

S = 1 - exp {-7 ( S) } (2.25) 

°o 1 


or 


4 C ~ O’ 

-4—-= exp { v( — ^) m } 

1- S a o 


.. (2.25) 


or log log 

1-S 


m log (a - a u ) - m log + log V 

.. (2.27) 


Eqn. 2.27 indicates that a plot of log log — (as 

1-S 

ordinate) Versus log (a - o^) as abscissa should be 
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a straight line for any particular volume Y - Y^, as 
shorn in Fig. 2.2. 


Further, the straight lines would he parallelly shifted 
with no distortion for different volumes of the same 
material, since a Q and m are independent of Y, In 
practice, the experimental data is first plotted as 
follows. The probability of failure S, corresponding 
to any stress level a can be found from 


S 


”0 

1 +1 


.. ( 2 . 28 ) 


where H = the total number of test specimens 

n a = the number of specimens that are ruptured 
upto the stress a 

While plotting the Weibull function is to be 

found by trial and error. First it is assumed that 

ct u is zero and if the plot happens to be a straight 

line, that implies the assumption a as zero is 

correct. If not a suitable value of o u is found by 

trial and error till the plot log log versus 

1-S 

log (0 - cr ) is reasonably a straight line. The 
slope of the straight line fomed earlier gives m, 
while the intercept on the abscissa gives c Q . Weibull 
applied his theory to a number of cases viz, strength 
of glass rods in tension, bending strength of porcelain, 
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tensile strength, of portland cement, filer strength of 
Indian cotton, breaking strength of cotton fabric strips, 
tensile strength of wood, stearic acid and plaster of 
paris, malleable iron castings etc. In all the various 
cases studied Wei bull verified the applicability of the 
linear fit of the distribution function with test data 
satisfactorily, L complete exposition of the various 
methods of fitting the experimental data suitable to 
^eibull's distribution and the related accuracies of the 
methods are given by Heavens and Murgtroyd (73). It is 
concluded by Heavens and Murgtroyd that information 
from small sets of data to a three parameter Wei bull 
distribution is not fully reliable j maximum likelihood 
estimates are found to be more accurate than lineari- 
zation and curve fitting. This study is of practical 
importance and is in agreement with the findings of 
Fruedenthal (14). 

2.3 FURTHER APPII CATIONS OF WEI BUhl'S THEORY 

Prompted by the work of ^eibull several investi- 
gations like Davidenkov et.al. (33), Weil and Daniels 

(74) , G-reene (46), Weil et.al. (47), Durelli and Parks 

(75) , and Evans and Pomeroy (53) examined the applicabi- 
lity of Weibull’s distribution function to several 
problems of static strength of materials. Davidenkov et^al f 



53 


studied the size effect on brittle fracture of steel in 
the light of w eibull 1 s theory and found the theory to 
be adequate in interpreting the size effect on mean 
strength and scatter. In addition the theory was also 
found to correlate the direct tensile and bending 
stresses at failure in that they are related by 

= (2m + 2) 1 / m 

- mean failure stress in bending 
= mean failure stress in tension 
= the material parameter 

yell and Daniels (74) derived the theoretical egressions 
for fracture probabilities in nonuniformly stressed 
specimens, such as beams under symmetrical four point 
loading, beams under centre point loading and beams 
under pure bending. The fracture probabilities in 
the ca.se of pure bending derived by the authors are 
verified by fitting by least squares technique the 
analytical distribution function to their test data on 
Columbia Resin specimens. G-reene (46) examined a possible 
generalization of'Weibull's distribution function in 
the following form 

N v 

S = 1 - exp { -V ( 2 c, a ) > .. (2.29) 

te =1 K 



CT t 


where a, 
b 

m 
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where c^. arc material constants. 

The distribution function with 

Z c k a k = 2.9 x 10’ 9 a 2 + 7.15 x 1CT 21 a 8 .. (2.30) 

was found to fit well with test data on pyrex Brand 
chemical glass specimens in bending, with one loading 
edge. As pointed by Greene, the expression relating 
mean strength to specimen size is extremely complicated 
even for the simple case of the series Z c^. where 
only two terms are present. Greene further -verified 
that tests on lime glass in bending with two loading 
edges are adequately characterized by Wsibull's distri- 
bution function. 

We il et.al (47) studied the applicability of the 
Weibull's theory to tests on Alumina. The applicability 
of the theory was satisfactory. The special feature 
of this particular study is the effect of parameters 
like surface treatment and t hemal history of the 
material on material parameters a u and m has been 
studied. It was noticed by the authors that grinding 
improves the value of m, leaving the value of ct u unaltered. 
Durelli and Parks (75) verified the satisfactory appli- 
cability of the linear plot of log (size) versus log 
(mean fracture stress), to the tests on Columbia resin 
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and plsxi glass in tension and tending. Richards (54,35) 
applied Weibullfs theory to characterize the scatter 
in the upper yield of mild steel in tension and tending 
and corroborated the applicability of Weibull's linear 
relationship between logarithm of stress at upper yield 
to logarithm of specimen volume ■ * Skinner (44) examined 
the applicability of both the theories of Weibull (48) 
and Frenkel and Kontorova (76) to study the scatter in 
testing of Anhydrite. Even though both the theories 
are equally applicable with in the range of sizes of 
specimens tested, the author remarked that Weibull’s theory 
is more suitable since the theory of Frenkel and 
Kontorova (76) is inapplicable at large sizes. This 
aspect will be further discussed in Sections 2.4 and 
2.8 of the thesis. 

The works of Evans and Pomeroy(53) and Salmassy 
et.al (49,50) on the application of Weibull’s theory 
to coals tested in compression and poreclain respectively 
are of special relevance to the thesis. The Weibull’s 
plots as given in the works of these authors are 
reproduced in Figs. 2.3 and 2.4jitcould be seen from 
the Weibull’s plots of cube tests on coal corresponding 
to various specimen sizes, that the plots are not 
parallel contradicting the assumptions of Weibull's 
theory. 'Also as remarked by Fruedenthal (14) a Weibull's 
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plot of the study by Salmassy et.al (49,50) on porcelain 
needs two different straight lines to characterize the 
strength distribution (refer to Pig. 2.5). The violation 
of the Weibull’s theory could be due to a change in 
fracture mechanism at different specimen sizes, other 
than that assumed in the basic development. This aspect 
will be further discussed in Section 2.8. 

2.4 THEORY OP FRENKEL AND KONT OROYA 

Using the weakest link concept Frenkel and 
Kontorova (76) investigated the scatter in strength of 
materials. The authors consider that the flaws are 
randomly distributed in the material. Further, assuming 
that the flaw strengths x, are considered to be statisti- 
cally normally distributed, so that 

i (x) = — exp < _ (x- v ) . } .. (2.31) 

/2i 3 2s 

where y = mean flaw strength 

s = standard deviation of flaw strength 

and § D (x) = probability density function of flaw strength 

Frenkel and Kontorova showed that the modal strength 
°m a SP 60 ^ 31611 °f size 7 containing n flaws per unit 
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volume Is approximately given "by 

a * y - s (2 log nV - 2 log 2 fit ) 1//2 ( 2 , 32 ) 

m 

where a = modal strength 

A more accurate formula for the modal strength 
a m can he derived (77, 78), using the results of 
Fisher end Tippet (79) and is given by 

a « v - s (2 log nV) l/2 + s 
m 2 (2 log nY) 1//2 

.. (2.35) 

In analytical procedure ofdetvtying the Eqn. 2.32 is 
discussed by Evans and Pomeroy (53). 

Apparently the Eqn. 2.32 of Frenkel and Kontorova 
(76) relating the specimen size and modal strength 
becomes invalid for very large or small values of Y. 

For very large values of Y the modal strength becomes 
negative, while for small values of V, due to the presence 
of log nV, the modal strength becomes imaginary. The 
reasons for the above discrepancies are that while 
deriving the result, the authors assume that n is very 
large, and the assumed normal distribution for flaw 
strengths associates a definite probability for negative 
flaw strengths. 
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2.5 THEORY OE El S HER AND HGDLOMON 

The statistical theory of fracture due to Eisher 
and Hollomon (80) is essentially based on the weakest 
link concept and considers the random flaw occurance in 
the material. The theory relates the strength distri- 
bution to flaw density and consequently to the size. 

The basic principle is that critical flaw in a popula- 
tion of flaws determines the fracture stress. Instead 
of directly building up a strength distribution function, 
Eisher and Hollomon start with an assumed distribution 
of flaw sizes given by the density function. 

p HH = e " e/hl •• (2-3+) 

n 1 

w here c = the crack size 

h^ = a constant, chosen crack length 

The function assumed given by Eqn. 2,34 indicates 
that the probability of flaw occurance decreases with 
flaw size. Further, it is assumed that the material 
fails according to Griffith’s (17) law so that 

s i* = A 1 c ~ 1//2 .. (2.35) 

where S. = the fracture stress 
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c 


A 


1 


Defining 


crack size 
a constant 



A iV 


• 1/2 


.. ( 2 . 36 ) 


Fisher and Hoi lorn on arrive at the flaw strength density 
function given by 

o -1/2 S ? 2 

P(S ? ) = — e ^ .. (2.37) 

^ 3 

S 2 

In order to take the random flaw orientation and number 
into effect, the authors consider that a specimen 
containing n flaws will fail when the normal stress 
on any one flaw reaches the critical value. The 
problem is studied by this concept for the co.se of a 
state of stress given by 


y 0 ; S-^. — Sy — or ; cc ( 1 *. (2.38) 

where X,Y,Z = there orthogonal directions 

S-£, Sy, S z = the. .srtrrooses in the three directions. 


The state of stress reduces to pure tension if a = o. 

Using the following notation 

* 


S 


r 



S Z 


P = 


s z 


. (2.39) 
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tc = / P(S ) dS 

oc x 


/ p(s ) {1 - 

a r M 


S - a 
r 

1-a 


} as. 


?(S)= — 

P 2 s 3 


i- s 2 

p 2 r 


.. ( 2 . 40 ) 

.. (2.41) 


wl lere 


S = the fracture stress of a specimen contain- 
n 

ing n flaws the authors show that 


P_ a* (?) 4P = n -i 11 " 1 1 — | ftp .. (2.42) 

n ’ a 5 13 

where P (P) = the probability that a specimen 

llj O- 

containing n flaws has a strength 
lying between P and P + dp . 

The corresponding relationship for a specimen under pure 
tension (<%*o) is 

p n,o* (p> = n < 1 ~ / p ( s r ) s r as r j 11 ’ 1 

o 

I f {P(SJ - /s } dS | .. (2.43) 

0 P o r r r 


For a given value of n , the above Eqn. 2,43 is to be 

evaluated numerically to evaluate P (P) , A plot of 

n y 0 

P n 0 versus P shows that the dispersion of strength 
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and mean strength decreases with increasing n. Assuming 

3 

surface flaws in glass occur with a density of 10 per 
sq.cm., the authors show good agreement between mean 
strength and specimen size. 

w hil e the theory of Fisher and Hollomon is conceptu 

ally elegant except for the a priori assumption of the 

distribution of flaw sizes, an evaluation of the fractupe 
probabilities numerically appears to be tedious. Also, 

If 

a correlation of specimen size with total number of 
flaws is not always easily possible especially when the 
flaws inside the material are effective in the fracture 
process as in the case of metals, as remarked by the 
authors themselves. 

2.6 THEORY OP KA3E 

Ease (81)in his study of the distribution of the 
tensile strength of rubber, considered that, of the very 
large number of flaws distributed throughout the material 
the most critical one determines the fracture stress. 

The author considers an equation of the foim 

S “ S o (1 *• .. (2.44) 

* 

where S » strength of the material in the presence of 
a flaw 

* 

^o = s ^ r ength of the material in the absence of 
a flaw 
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a* = a constant 

x. = flaw size 
1 

With the above Eqn. 2. 44, the problem of the study 
of the distribution of strength S reduces to that of 
the distribution of the flaw size . Ease assumes 
that the flaws are distributed exponentially, i.e. 

f ( Xl ) = X e " .. (2.45) 

where f(x^) - the probability density function of 

flaw size, 
i = a constant. 

Corresponding to this density function the most 

* 

probable maximum flaw size x in a specimen of volume 
V is shown to be 

* - l.Qg Vn 
x ~ X 

where n * the average number of flaws per unit volume. 
Consequently the modal strength of a specimen 

** 

S , having a volume V is given by 

S «S t (1.« losf!) .. (2.46) 

Assuming that the number of flaws n is very large, Ease 
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further shows that the strength distribution is given by 
E (n ) = exp {- exp (-n ) } .» (2.47) 


where 


x 

11 = -TTTf 

S a 
o 


(s - s ) 


.. (2.48) 


E (n) = strength distribution function, 


X *)f 

So that, knowing —jr* and 5 one can readily get the 

Sa 

distribution of strength S . Experimental investigations 
yield the values of mean strength <T and standard 
deviation S^. Corresponding to the above distribution 


function it has been shown that 

S** -- 

‘ X 


-if * 


. 0.1577 


.. (2.49) 


* * 

S’ a 

o rc_ 

X 1(6 


= S 


d 


.. ( 2 . 50 ) 


"7ith a knowledge of c and S^, the function E(n) in 
Eqn. 2.47 can be characterized. Ease applied his theory 
to verify the scatter in his tension test results on 
rubber and found the predicted form of E(n) to be 
agreeable with the experimental one. 


2.7 OTHER PHEHCMMQLOGICAL THEORIES 

The statistical theory of strength due to Valkov 
(82) is a development of the basic concepts of continuum 
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mechanics to random micro-media or a microscopically 
inhomogeneous medium. In contrast to the theories 
reviewed earlier, the theory of Valkov is deVoted to the 
determination of stress fields in a random medium and 
therefore is a study of statistical continuum theory 
analogous to the work of Beren (102). Brady (83) advanced 
a failure criterion for rock like materials in which 
microcracks develop and propagate under increased 
loading. Considering the microcrack density at various 
stages of loading upto failure and the rupture of 
material according to Mohr - Coulomb criterion, Brady 
advanced a failure criterion that, fracture occurs when 
the total volumetric strain due to microcracks attains 

a critical value. The theory is shown to be in quali- 
tative agreement with various aspects of mechanical 
behaviour of rock like materials and is used to explain 
the nonlinear stress strain curve as affected by 
micros racking. 

It may be noted that the two theories mentioned 
in this section, eventhough are termed as statistical 
strength theories, treat altogether a different problem 
from that treated by Teibull and others. The team 
’’statistical aspects of strength and fracture" is used 
for the study of creep - rupture of concrete as a 
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stochastic process "by Hori (129 ), The aforesaid theories 
do not deal with strength specifically as related to 
size although the mechanical behaviour is explained on 
the basis of random internal structure of the material. 

2.8 DISCUSS I Oh 

The application of weakest link concept in the 
study of statistical strength theories of materials 
is popular because of the simplicity in mathematical 
formulation of the problem. However, the implications 
of the weakest link concept are vulnerable to some 
criticism* In the study of a continuous material in the 
presence of randomly distributed flaws, the weakest link 
concept centres attention only on the most critical or 
potential flaw and the interaction between the various 
flaws is neglected. Qualitatively, the essence of the con- 
cept is that larger the specimen size? greater is the 
probability of a critical flaw occuring. As such, a 
majority of the studies like that of Weibull (48) and 
Frenkel and Kontorova (76) have been shown by Epstein 
(78,84), Saibel (85) and Romauldi (86) to belong to a 
class of problems in exitreme value statistics, in which 
the distribution of the smallest value in a sample taken 
from a population is studied. Apart from the statistical 
distribution of flaws in the material, the various theories 
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based on weakest link concept do not take care of any 
of the aspects of actual mechanisms of fracture in all 
types of materials.. Table 2.1 gives a, comparative 
tabulation of the assumptions in the mathematical 
formulation and the final results of the various 
theories discussed in the preceding sections. 

2,8.1 Comparisons and Limitations of the above Theories 
A general limitation of the various theories based 
on weakest link concept reviewed above is an adhoc and a 
priori assumption of the form of strength distribution 
function, flaw strength distribution or flaw size 
distribution. Eventhough an assumption of some kind 
is justifyable on the grounds of simplicity, the real 
problem of identifying a distribution to fit the 
scatter in materials strength results obtained by actual 
tests is circumvented. Consequently, the true fracture 
mechanisms are wiped out and therefore it should at best 
be considered as a hypothesis or an assumption which 
should be checked while trying to applyfra particular 
material. By such assumptions, the characteristic 
features of the scatter of the material are a priori 
fixed by the assumption made. 

Inspite of its remarkable simplicity the distri- 
bution function assumed by w e ibull (48) automatically 
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fixes the form of size-mean strength relation. It can 
he easily observed that the actual mechanism at fracture 
can he different at different specimen sizes, like 
for example, depending on the energy content for 
driving a potential crack. If the specimen is large, 
the total energy of the material under load is higher 
and consequently .the crack driving force is more, than 
that of a smaller specimen. The larger specimens with 
a potential crack are likely to fail in a brittle 
manner as compared to smaller ones. Consequently, it 
may not he always possible to characterize the size 
mean strength relation by a simple a priori fixed form. 

Frenkel and Kontorova (76), instead of assuming 
a distribution function directly assume that the flaw 
strengths in the material are distributed normally 
and that the density of flaws is very large. As a 
consequence the size-modal strength relation obtained 
by this theory is inapplicable at very large or small 
volumes. Physically it can be seen that at very large 
sizes the strength should approach a finite value and 
also at smaller volumes the theoretical strength of the 
material in the absence of any flaws. 

Fisher and Holl onion (80) assume flam length 
distribution in the material and relate scatter in 
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strength with the density of flaws, Eventhough the 
results of the theory are not available in closed form, 
it can he expected that the assumed form of flaw 
distribution influences the form of the strength 
distribution as well as size-mean strength relationship, 
The theory of Ease (81), assumes the nature of flaw 
distribution as well as the relation between strength 
and flaw size, ihe applicability of the relationship 
between modal strength and specimen size of Ease is 
limited at very large volumes. 

Of the many theories referred to above, that of 
teibull (48) is applied to a wide variety of problems^ 
However, as could be seen the applicability of the 
theory is limited in the case of some materials which 
will be discussed further in the thesis, 

2,8.2 Scope for Alternative Approach 

Inspite of the assumptions implicit in the weakest 
link concept viz, that only the potential flaw deter- 
mines material strength and flaw interaction effects are 
neglected, the application of the concept is still 
welcome because of the mathematical simplicity in 
studying the problem of characterization of scatter in 
materials testing. However, as could be noticed from 
the theories reviewed earlier any a priori assumption 
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ox the fom of strength distribution function or flaw 
size distribution essentially limits the scope, in 
that the form of size mean strength relation is fixed. 

As could be seen from the studies of Evans and Pomeroy 
(53) on coal in compression, the Peibull plots at 
various specimen sizes, Pigs. 2.3 and 2.4,- event hough 
straight lines, are not parallel. This implies that 
the parameter r m’ in the lei bull ' s distribution that 
determines the slope of the straight line in the leibull’s 
plot is a function of specimen volume. In effect this 
means that the assumed form of distribution functien 
is incorrect. Similarly the results of Salmassy St.al 

(49 j 50) on porcelain indicate the need for two straight 
lines in leibull's plot, Pig.2.5,to characterize the 
scatter as also remarked by Fruedenthal (14). These 
remarks show the limitations of applicability of the 
above mentioned formulations for all solid materials 
eqrelly well, because of the inherent assumptions, A 
possible approach to by pass these problems is to seek 
a methodology by which a strength distribution function 
can be arrived at formally, suitable to a particular 
material as given by its size-mean strength relation* 

In this manner, the true material characteristics can be 
reflected in the formulation* A formulation of this 
problem and its solution in special case is discussed in 
the following chapters. 
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CHAPTER THREE 

PROPOSED POHMULATION OP THE PROHklM OP 
~ "S‘lZE~*EPPEO T OH ST'REN'GTH ~OP MATElgSE5 “ 

3.1 BASIC MOTIVATION POR THE PRESENT P OMUL ATI ON 

Prom the review of various statistical strength 
theories in the earlier chapter, one can notice that 
the weakest link concept, even though an idealization of 
material's mechanism of fracture, does provide a basis 
for mathematical formulation and application. Further- 
more, the various theories discussed in Chapter 2 start 
with an a priori assumption of the form of strength 
distribution function. As a consequence of such assump- 
tion, the form of size-mean strength relation or size- 
modal strength relation is automatically fixed. The 
accuracy in the application of any particular theory 
lies in the problem of how best the material under 
consideration can be characterized by the size-mean 
strength relation or size-modal strength relation 
corresponding to the applied theory. However, as 
mentioned in Section 2.8.2, application of Weibull's 
hypothesis to data available caCcaL and porcelain showed 
the limitations of the basic formulation in taking 
care of behaviour of these materials., A plausible way 
of finding a more suitable distribution function 
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reflecting observed material behaviour, could be made 
by constructing the same from first principles taking 
the size-mean strength relation as the characteristic 
of the material and getting the distribution function 
corresponding to the material rather than assuming the 
form of the distribution function a priori. In the 
present Chapter, possibilities of this approach are 
examined in detail. 

The motivation for the present work is that of 
Tsai and Kolsky (87) in which they examined the apparent 
fracture strength of glass plates as dependent on the 
diameter of the indentor. The authors constructed a 
theoretical strength distribution function corresponding 
to the mean fracture strength - indentor radius relation 
as given by experimental investigations. The authors 
concluded that scatter prediction corresponding to 
large indentor strengths is satisfactory by the statistical 
approach. 

In this chapter, the close analogy of the problem 
of size effect on strength and that esflndentor size on 
fracture strength of glass as attempted by Tsai and 
Kolsky (87) is used for formulation of size effects 
on strength of materials which is an alternative 
approach to presently available theories (as pointed out 
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in Section 2.8). 

3.1,1 Spec ifi cation of th e P roblem 

In the following, the problem of constructing 
the strength distribution function that corresponds to 
a given form of size-mean strength relation (obtained 
experimentally) is attempted. The distribution function 
is required to reveal the general features of mechanical 
behaviour from typical tests on materials for strength 
with size as a parameter, viz, the decrease in scatter 
with increase in specimen size. The method of approach 
and solution is developed for this problem in a manner 
similar to that of Tsai and Kolsky (87) in their investi- 
gation of tensile strength of glass as dependent on 
indentor diameter. 

3.2 ASSUMPTIONS MADE IN THE FOMULATION 

The basic assumptions of the nature of failure of 
material are the same as those of yeibull, viz, that 
the material obeys weakest link concept in that the 
fracture of any typical element is the same as that of the 
total material; the flaws initiating rupture are 
randomly distributed throughout the material and the 
interaction between various elements with flaws is 
neglected. The most potential flaw that forms the 
weakest link determines the rupture of material. 



75 


3.5 THEORETICAL MODEL OF THE PROBLM OF SIZE EFFECTS 

OH STRENGTH 

As has been discussed in Section 2.1 of the thesis, 
a general feature of the tests on strength of materials 
is the decrease in mean strength and scatter with 
increasing size. The size could correspond to volume, 
area or length depending on the nature of material and 
their testing. In the following presentation, the 
volume of the material is considered to represent size 
and a change to area or length, if need be, depending 
on the situation does not pose any problem. 

Considering the strength of a specimen of size V, 
let the total volume be divided into several infinitismal 
elements of size AY, Let the probability, that a 
potential flaw exists in an elemental volume AY, and 
the flaw is of such size that fracture occurs when the 
applied stress is between a and a + da , be defined 
as {<}> (a ) dc> a V. It may be noted that this 
definition is analogous to that of the occurance of an 
event in an interval At in a Poisson process of 
parameter X which is given by Xdt. In the present 
problem, it so happens that the corresponding parameter 
is dependent on the applied stress a and is given by 
<f> (a) da . Having defined {<{> ( a) da } A Y as above, 
it can be written that the probability P(o) that 
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fracture occurs in the element under consideration for 
any stress level lying between 0 and a he given by 

P(<0 = { f° <t>(cr) da > x AT . . ( 3. 1 ) 

o 

In many materials the size-mean strength relation 
converges asymptotically to a value , the lowest 
mean strength and this can be incorporated in the 
analysis by specifying 

$(<0*0 for a < a L .. (3.2) 

_ *, . 

Using Eqns. 3.1 and 3.2 the probability P (a) that 
fracture cannot occur in the elemental volume AV for 
a stress level lying between 0 and a is given by 

P*(<0 = 1 - { f° 4> (a) da > AT = 1 - G(o) AV 

°Jj 

where G- (a ) « f a $ (a ) da .. (3.3) 

a L 

If the total volume V of the specimen is considered to 
be divided into several infinitismal volumes A"V\ 

{ i = 1,2 ... co} f the probability that failure 
cannot occur anywhere In the total specimen $ (o) 
is given by 

$ (a) » nil- G(o) AT i j ; 1=1,2, 3 ... 

« * ( 3.4) 
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It may "be noticed that in writing Eqn. 3.4, the weakest 
link concept is implicitly used in the sense that the 
failure of total material Is considered to he same as 
that of any one of the elements A failing. 

The infinite product in Eqn. 3.4 may he simpli- 
fied using the following inequalities. Eor any x^, 
so that 0 < x^ < 1 one can write 

nr 

e" x i > (I-Xjj > e 1 - 1/2 .. (3.5) 

-x. -2 x. 

or xO-om) < % e 1 » e 1 . . 


jc • x 

and > e 1 it( 1-1/2 x. 2 e i ) 


.. (3.6) 


In the limit, as x -* 0, the above inequalities Indicate 


that 


n ( 1 -x^ ) - e i 


.. ( 3.7) 


Identifying x i as G(or) A y. i n the Eqn. 3.4, one can 
write 

, -G(cr) Z A V. 

% (I-G(ct) AV i ) = e 1 


as a Vi -*• 0; 
n (l-Gr(ff) A V ) 


-G(cr) / dV -G(ff) 7 
e V = e 

-V f +(a) da from Eqn. 3. 3 
e °L 


* « 


(3.8) 



78 


Since $(<r) is the probability that failure does not 
occur anywhere in the specimen of volume Y, the proba- 

■¥r . 

bility $ (a) that fracture occurs for any stress 
lying between 0 and a , is given by 

4%) = 1- * (a) * 1- e~ V 


-Y f° <K CT ) dc? 

- 1- e a j, for a > 

.. (3.9) 

= 0 for cr < cr 


The density function f( a ), corresponding to 4 (cr), 
namely the probability that fraeture occurs for an 
applied stress level lying between a and a + do is 
given by 


f (O 


d (cr) 

d a 


.. (3.10) 


Also, if a is the mean fracture stress, one has 

oo 

_ £- a f <°) a 0 » ' 

a ~ ~ •/<tf(o)do»/af(«)dff .. (3. 11 ) 

/ f(o) do 0 °L 

o 

Since f( a ) * 0 for a < a _ 

— JLi 

Integrating Bqn. 3.11 by p aits and evaluating the limits 
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a 


CO 




j_ a $* (cr) - / $*(<?) dcr j 


oo 


= { 0 $ (cr) } 

0 

L 




<J>(0 d5 „ 

} da 


or 

CO 

7= a + / 
h 0 

L 



<KO dC > 


da 


where £ is a dummy variable of integration. 


( 3 . 12 ) 


The above Eqn. 3.12 is the relationship that 
relates specimen size V to the mean strength o’. If 
<f> (O or 'K 0 ) is evaluated corresponding to a given 
size-mean strength relation, one can find f(a ) which 
is the density function denoting the probability of 
failure, of specimen for a stress lying between a and 

-ft 

a+da and $ (a ) which is cumulative probability of 
failure using Eqns. 3.9 and 3.10. A method of solution 
of Eqn. 3.12 will be discussed in detail in the next 
Chapter. Itr may be noted that Eqn. 3.12 is a non- 


linear integral equation 
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3.4 DISCUSSION AND COMPARISON WITH THE EXISTING- THEORIES 
In contrast with the methods of any of the 
statistical strength theories (Chapter 2), the present 
approach seeks a method by which the strength distri- 
bution function can be constructed corresponding to a 
given form of size-mean strength relationship, with 
no a priori assumption of the form of strength distri- 
bution or flaw distribution function. The present method 
in effect is converse in its method of approach to exist- 
ing ones, in that the distribution function is obtained 
subsequent to a given size-mean strength relationship 
while in all the theories of Weibull. and others, some 
assumption of strength distribution is made (Chapter 2) 
and size-mean strength or size - modal strength relation 
is obtained which is fitted with the experimental data* 

While the ease and simplicity with which the 
theories of Weibull and others can be applied to some 
materials, are undispu table, it is imperative that alter- 
native methods of approach should be available to treat 
cases of materials that do not suit the available theories 
like for example the coal cubes and porcelain results 
discussed in Section 2*8. 1 r wherein the Weibull plots 
corresponding to various sizes are not parallel for 
coal cubes and two straight lines are required to fit 
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data <£n porcelain. While a possible approach is to try 
alternate forms of distribution functions, it may be 
relatively easier and element to try alternate forms 
of size-mean strength relationships corresponding to 
the actual experimental data and the problem formulated 
is aimed at dfcvieing a general method by which distri- 
bution functions can be obtained, given the size-mean 
strength relationship. 



S3 


a 


{ 1 + 



.. (4.3) 


where a ^ is the mean lowest strength; K and n are 

material parameters. 

Eqn. 4.3 chosen above indicates that as 7 increases 
o’ decreases and tends to a asymptotically as 7 
tends to he very large. 

Substituting Eqn. 4.3 in 4.1 


a L K 



°° -7 f° <f»(0 d? 

/ { e > do 


.. (4.4) 


To simplify the above Eqn. 4.4, let the following 
notation be used 

/ <fr(S) dS = Ku) .. (4.5) 

a l 

so that 

a L K 


Differentiating Eqn. 4.5 under the integral sign 

= ( u ' = *00 .. ( 4 . 7 ) 

du 


/ 

a 


-7 4(u) 


du 


(4.6) 


so that, knowing i|>(u), <J>(u) can be evaluated by 

differentiating if>(u). 
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To evaluate ^(u), let 


l>(u) 

and u 


t 

u (t) 


so that du = u f dt 


.. ( 4*8) 

.. (4.9) 

.. ( 4 . 10 ) 


Assume that as u t — a o and as u — t — °® 

where a 0 is any constant greater than or equal to zero. 

^ith this notation, Eqn. 4.6 takes the form 


a L K 



r -Vt . 

J e u ' dt 

a 


.. (4.11) 


or equivalently, hy dividing either side hy V 


„ - u' dt * a T K V" (n+1) 

a n V L 


/ Q 

: o 


.. (4.12) 


Taking the Inverse Laplace Transform of Eqn. 4.12 with 
respect to V one obtains 


/ u* H(s-t) dt 
a. 


c^K s 


n 


.. (4.13) 


-o r ( n+ 1 ) 

where H( s-t) is the Heaviside unit function given by 

H(s-t) = 0 for s < t 

* 1 for s > t 


.. (4.14) 


s ~ is the parameter in the Laplace Transform 
and r (n+1 )« gamma function of argument (n+1) 
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Using the properties of H(s-t) given by Eqn. 4.14, 


Eqn. 4.13 takes the foam 



T r n 

s a x m s 

f u ! at = — — — 
a 0 r (n+1) 

.. (4.15) 

or 

o t K s n 

u(s) - u(o^* 

r (n+1) 

.. (4.16) 

Eqn. 

4.16 can be revcritten as 



•t Tf + n 

u(t) = u(c^) + — 

* r ( n+ 1 ) 



. K t n , 

— 0 T { 1 + 11 1 } 

L r ( n+ 1 ) 

.. (4.17) 


It may be verified from Eqn. 4.17 that the assumption 
made in introducing u, viz, that as t — u -*• 00 is 
satisfied. Erom Eqn. 4.8 

4 (u) = t 

Eqn. 4.17 can be rewritten as 


u , {1+ , 

L r ( n+ 1 ) 


... (4.18) 


so that 

♦fa) > 1/n 

K a T 


.. (4.19) 


Since from Eqn. 4.7 
0 (u.) - 4’ (u) 
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$ (u) can "be evaluated using Eqn. 4*19 as follows 

r( .4N 1/a d/n}-1 

+ (a) = (u) } . _ . ( u _ a L ) (4.20) 

K a L 


From Eqn. 4.20 and the restriction on *(o) for a < 
(See Eqn. 3.2 of Section 3.3 in the previous Chapter) 


one has 


♦ to = . i . ( 

T r „ 41 

1*. 0 -r 


O- Oj.) 


(l/n)~1 


* 0 


G > 0. 


a i °L 


.. (4.21) 


Eqn. 4.21 is the solution of Eqn. 4.1 corresponding 
to the size-mean strength relation of Eqn. 4.3. By direct 
substitution it can be verified that <K a ) from 
Eqn. 4.21 satisfies Eqn, 4.1. The expressions for 
$ ( CT ) and f(tf ) are given below using Eqns. 4.21, 3.9 
and 3.10 of the preceding chapter. 

1/n * / 

(o) = i_ exp [ -y { j 


= 0 


for o > cr j 
for a < cft 


(4.22) 


1 (a) . I {_£Istl. j . (a-^p/n-l 


esp l. T <I£aUl > 1/n . (o _. ) 1/ “ ] fo r0 >. : 


( 4 . 23 ) 


f or CT 1 % 
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Corresponding to $*(<*) and f(cr) given by Eqnr/. 4.22 and 
4. 23 above, the standard deviation is given by 


S = { / G^f ( a ) da 

a o 


# ) 1/2 


Si GTEsII . 1 

T n i rpn+1) 


where S, = standard deviation 
d 

The coefficient of variation C is given by 


. (4.24) 


C 

T 


o + f) 


n 


I (2 n+1) _ 1 

rVD 


. (4.25) 


-.1,1 Speci alization to itei bu ll's Resul b 

^eibull's strength distribution function is given by 


. -V (~“) m 

4 (? ) = 1 - e o 


.. ( 4 . 26 ) 


and the size-mean strength relation and the variance S^ 
corresponding to this distribution function are given by 


a — a + 
u 


Y 


1/m 


r ( 1+1/m) 


.. (4.27) 


= a 2 Y~2/m { r (i+2/m)- r 2 ( 1 +i/m) } (4.28) 


The size-mean strength relation considered in the present 
formulation is given by 


° = a T {1+ 

1 yR 


.. (4.29) 
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and the variance is given by the square of standard 
deviation o o 


(S d )' 


8 L g 

V Zn r^(n+l) • 


_ D 


.. (4,30) 


Comparing Sons. 4.29 and 4.27 the following correspondence 
between parameters in Ueibull’s equations and those in the 
present formulation can be established 


a a 


K ^ c ^(1+l/ra) and n -v i/m .. (a. 31) 


1 u ’ 1 " o 

Using the relations 4.31 in Sqn. 4.30 

li Il+i/jai { . r( 

Tr2/irL r^U-l/m 


(S a )' 


- 11 


— a 


r £ -(i+i/m) 

' y*"2/m { r (i+2/m)- r 2 ( l+l/ m ) > .. (4.32) 

which is the same as ^ei bull’s result 


4.1.2 Speciali z ation to Bolotin's Result 

In the notation of Bolotin (72) the coefficient of 
variation in strength as affected by specimen size is 
given by (Refer Section 2.2.2) 

b /-) 1/a ♦(«) 


c 

v 


V. 


.. (4.33) 


, , f On 1 /a 

a + b (y) 


where <J> (oc)= 


V 


- r m±.?Z“U. _h 

O l 


.. (4.34) 


r ( 1*M /oc) 

The size-mean strength relation of Bolotin is of the form 


T 


* = \ { a + b (^2) 


On 1/a 


V 


} 


.. (4.35) 
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where 


R 


R 0 ,a,h,a 

T_ 


the mean strength 
are constants 
some standard volume 


The corresponding results of the present approach are 


v 


( 1 - 


V 


n 




l ( 2nt 1 ) _ 

r 2 (n+1) 


.* (4*56) 


and 5 = 


a l f1+ ^ } 


•. C4*3T) 


Comparing Eqns. 4.35 and 4.37 the following correspond- 
ence between various constants can he noted 


a R, 



a h -*■ R h V 1//a 
L oo 


and n i/a 


The coefficient of variation given by Eqn. 4*36 > with 
the above changes in notation is given by 



V*l/ a D 


1 + 


» Rq * cl 




1 7a 


o 


r O+lM 

r (1+1/a) 


.. (4.38) 


b(Y 0 ) 


i/a 


u 7 0 1 /“ + y 1 /“ 


J 




Os 1 /a 


T 


a+b(^) 


os 1 /a 


a 



( 1+2 /a) 
2 ( 1+l/a) 


r ( 1+2/a) 
r 2 ( 1+1/a) 


-1 


1 


Which is the same as that given by Eqn. 4*35. 
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4.2 FEASIBILITY OF CLOSED FORM SOLUTIONS 

The solution of Eqn, 4.1 for <K°) corresponding 
to a given form of size-mean strength relationship will 
he in closed form whenever the Inverse Laplace transform 
in Eqn. 4.12 as well as an explicit -form of ^(u) in 
Eqn. 4,18 can he obtained. However the above restri- 
ctions are stringent in that both may not be obtained 
without resort to numerical procedures. Tsai and 
Kolslcy (87), in their attempt to get a fracture stress 
distribution function (corresponding to a given indentor 
radius-mean fracture load relation), obtained the distri- 
bution function in closed form only for a quadratic 
relation between indentor radius and fracture load. 

In the following, the solution of Eqn. 4.1 for 
<K a ) corresponding to a more general form of size-mean 
strength relationship given by 

a = d (1+ MU } (4,39) 

D(V) 

where D(Y) is a polynomial of degree n 

N(V) is a polynomial of degree n-1 or less 

is attempted. The form of Eqn. 4.30 reflects that as 
V increases o decreases and approaches cr^ in the 
limit. The following procedure is analogous to what 
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has "been followed in the earlier section. 
Substituting Eqn. 4.39 in Eqn, 4.1 


° r IT(Y) °o { -Y f° 4(£) d£ } da 

/ .e a -r 

L(Y) a 


C L y(T) 

D(V) 


00 -Y 4(u) 

/ e du 

c 

1 


and 


a L h(V) 

YD( V ) 


/ 

cc 

o 



u' dt 


.. (4.40) 


. . (4.41') 


.. (44?-) 


where 4(u), t, and oc q are defined exactly as before- 

Taking the Inverse Laplace transform of Eqn. 4.4%. on 
either side with respect to Y 

«> . -1 a¥(Y) 

f u’ H(s-t) dt * L {— > .. (4.43) 

a YD(Y) 


s -1 a T lT(Y) 

/ u f dt = L {— - } (4.44) 

a YD(Y) 


— 1 

where L denotes Inverse Laplace transform. The 
quantity on the right hand side of Eqn. 4 .44 can be 
evaluated as (68) 


-1 


L 


a h(Y) 

{ 

YL(Y) 


} 


n+1 

l 

xo=1 


2 > } 

35*^7 



.. (4. :-5) 
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where a,, a 2 , ... a n+ . 1 are the roots of the equation 


ai 


prom 


u(s) - 
or 

u(s) 


T> (Y) 

= Y D(V) 

=* 0 



i 

I) denotes the d 

erivative 

of 

the polinomial 

5 . 

Eqns . 

4. d4 and 

4.45 





» u( s) - 


n+1 

v 

tfCam) a m s 


- u(cQ 

0 L = a L 

l 

m=1 

i 0 < 

. * (4.46) 


n+1 

{ 1+ I 

f 

cl 

in 

e 

s 

} 

(4.47) 


m=1 

15 





or 


U- 0 T 
{ !■} 


T \ ♦<»> 

\ TO e 

m=1 


. . (4.4 8 ) 


Eqn. 4.48 above gives the functional form of ^ (u) which 
when evaluated could be differentiated to get 4>(c). 

It can be seen that in a general case, 4>(u) cannot be 
explicitly written as a function of u in which case 
a numerical approach is essential in proceeding further. 

4.3 MERE CAL TREATMENT IN MORE GENERAL CASES: SUGGESTIONS 

In the example considered in Seetion 4.2 it was 
noticed that 4(u) that is to be evaluated (which when 
differentiated yields <j>(u) or 4 >(ct ) , could not be 
obtained explicitly. In such cases, the relation 
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between ifi(u) and u has to bo numerically tabulated using 
Eqn. 4.48. The tabulated function when differentiated 
would yield <!>(u) and further evaluation of $ (° ) 
and f(a) according to Eons. 3.9 and 3.10 of the preceding 
chapter is also to be carried, out numerically. A resort 
to such numerical evaluation procedures should be 
commensurate with the accuracy warranted by the nature 
of the problem under consideration as compared to avail- 
able net nods like that for example that of Weibull. 

4.4 DISCUSSION 

The present formulation and method of solution of 
the problem of size effects on strength is an attempt to 
pave the way to construct distribution functions suitable 
to a material behaviour as given by its size-mean strength 
relationship, The feasibility of such an approach is 
shown by getting the distribution function of "Weibull 
starting from the corresponding size-mean strength relation- 
ship in a reparameterized form. However, construction of 
strength distribution functions suitable to more general 
forms of size-mean strength relations needs numerical 
approaches as was exemplified by the problem considered 
in Section 4.2. The scope of the present thesis is 
only to examine the feasibility of obtaining strength 
distribution function starting ’with a given form of 
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size-mean relationship with weakest link concept as 
basis, as this approach is more general in its appli- 
cation. Consequently no attempt is made to solve numeri- 
cally the problem of size effects on strength for more 
general formso^size-mean strength relationships. 

Further investigations are necessary to obtain the 
strength distribution functions both in closed form as - 
well as numerically, the size-mean strength relation- 
ships being obtained from tabulated data on strength 
of materials with specimen size as the variable. 



CHAPTER El YE 


CHARACTERIZATION OE DIRECT TENSILE STRENGTH OF 
' 0 'ONffREfS~BY - A- ' PI 5TRE BU TI ON FUNCTION 


5.1 INTRODUCTION 

Concrete is a very important engineering material 
"but an understanding of its properties and mechanical 
behaviour is still inadequate. Employing the classical 
concepts of continuum mechanics very little can be 
achieved in understanding the large scatter associated 
with concrete or the peculiarities in its mechanical 
behaviour. Considering the complex internal structure 
of concrete, it is easy to visualize the inherent 
difficulty in arriving at a suitable formal mathematical 
model for concrete. Inherently, concrete is a hetero- 
geneous material consisting of mortar and aggregate 
phases mainly. The interface between the two phases 
being a zone of relatively weak bond, microcracks 
are present even at no load (89), due to shrinkage. 

The microcracks or interface bond failures act as stress 
raisers and cause progressive failure of the material 
under loading. While it is certain that concrete can . 
hardly be idealized as brittle material because of 
the various forms (21), other than surface energy, 

r 

in which the energy can be dissipated during crack 
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propagation, quantification of the various energy 
dissipation mechanisms appears to he hopcle-ssly difficult 
at present. 

In the present chapter an attempt is made to 
characterize the size dependent strength and scatter of 
concrete in direct tension, hy a distribution function 
using the results of the previous chapter. Eventhough, 
concrete cannot he considered as a material that closely 
obeys weakest link concept because of the presence of 
various energy dissipating mechanisms, the decrease in 
scatter and mean strength with size prompts one to 
examine the applicability of the distribution function 
obtained based on weakest link concept. 

In the following, a brief outline of the works of 
various investigations that brought to light the flaw 
sensitivity in concrete and the applicability of the 
concepts of fracture mechanics to concrete is given, so 
as to provide a basis for the subsequent development in 
the chapter wherein a study of the strength distribution 
of concrete is made. 

5.2 BEHAVIOUR AND STRENGTH OE CONCRETE 


The earliest proponent of the application of 
Griffith's theory to concrete is Kaplan (90), whose 



worlc gave qualitative evidence of the applicability. 
Following the work of Kaplan, GluckLich (91), explained the 
difference in compressive and tensile strength of concrete, 
higher strength associated with smaller aggregate mixes 
and other related phenomena using the concepts of 
fracture mechanics. G-luckLich explained the severity 
of the presence of microcracking around crack tips in 
concrete as well as the effect of aggregate phase in 
influencing the energy demand and dissipation curves during 
crack propagation in concrete. Romauldi and Batson (92), 
Sridhar Rao and Parimi (93) established the role of 
fibers in arresting crack propagation in fiber reinforced 
concrete. The works of the above investigators as well 
as those of Lott and Kesler (94), Welch and Kaisman (95) 
have established that fracture mechanics concepts can be 
used to study the material properties of concrete. 

Further evidence regarding the effect of microcracking 
in concrete on the fatigue behaviour, cumulative damage 
and sustained load behaviour is due to Raju (96), 

Neal (97), Sushil Chandra (98). 

Because the microcracks in concrete are random 
in their distribution in the material, scatter in test 
results is often present. However as has already been 
remarked, concrete cannot be rationally idealized as a 
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"brittle material since catastrophic crack propagation 
and failure is not truely present as in materials like 
glass. Tension failures in concrete are relatively 
sudden and spontaneous as compared to compression tests. 
This is because the tensile stress field is more 
favourable for an arbitrarily oriented flaw to propagate 
as compared to compressive stress field, on the basis 
of which Glucklich (91), Interpreted the ratio of the 
compressive to tensile strength of concrete. Further 
an examination of typical tests on concrete for tension 
(45) and compression (41) indicate that the mean strength 
and scatter decrease with increase in specimen size 
displaying that the features are akin to a system 
obeying weakest link concept. Applicability of 7/eibull ’ s 
strength theory to concrete has been discussed by 
Popovics (99), who suggested that strength ratio of 
two specimens is related approximately to their volume 
ratio as given by 




m 


. . ( 5 . 1 ) 


where cr ^ = mean strength of specimen of size 

o~2 = mean strength of specimen of size Y 2 
and m = a constant, the material parameter. 
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It can however he noticed that such a strength ratio 
as given by Eqn. 5.1 implies that the rate of decrease 
of strength with size is independent of specimen volume 
as well as an indefinite strength decrease with increase 
in size. This is a consequence of fact that the lowest 
mean strength of the material which is asymptotically 
approached as the specimen size increases has been 
assumed to he zero. Bolotin (72) discussed the appli- 
cability of the three parameter strength distribution 
function as given by Eqn. 2.16 discussed in Section 2,2.2 
of the thesis, to the compressive and flexural strengths 
of concrete. Such a. three parametered distribution 
function is essential to accomodate for the existence 
of non -zero lowest mean strength of a material. 
Consequently an attempt is made to characterize the 
scatter in direct tensile strength of concrete, Tby a 
distribution function using the results of the previous 
chapter, 

5.3 E7ALUA1I GST OP MATERIAL PARMBfERS 

In the following an attempt is made to find the 
material parameters , K andn as applicable to 
direct tension results of concrete, in the equation 

o' — o (1+ ~r } 

L yll 


• • 


(5.2) 
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where c = the mean strength 

Y - specimen size 

A comprehensive study of the effect of specimen size 
on the direct tensile strength of concrete is due to 
Kadlecek andSpetla (43) and it is proposed to use the 
results of the authors in the present investigation. 

/ V 

Kafllecek and Spetla have eliminated the effect of 
slenderness ratio, viz the height to lateral dimension 
ratio and the results provide the volume effect on 
strength as observed in direct tension testing. 

5.3.1 Details of Ex p erimental Data 

Kadi ec ek and Spetla (43) have conducted fasts in 
direct tension on cylinders and prisms; two series in 
each case with specimen volume as the variable. The 
details of the test specimens and the experimental results 
are given in Table 5. 1 and Digs. 5. 1 and 5. 2. The last 
column of Table 5. 1 contains the size-mean strength 
relations constructed in the present investigation and 
are explained in the following section. The Pigs. 5.1 
and 5.2 depicting the test data are reproduced fr~>m the 

v V 

paper of Kadlecek and Spetla (43). 

5.3.2 Material Parameters 

V V 

Kadlecek and Spetla suggested the size mean 
strength relations as given in Table 5.1 to be 
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T ABL E 5.1: DETAILS 0E EXP ERIMENTAL DATA OE DIRECT TENSION 
TEST RESULTS ON CONCRETE 


Size-Mean Strength Relation 


Test Series 


Originally giver 
by gadlecek and 
Spetla (43) 


Constructed 
in the pre- 
sent Investi- 
gation 


CYLINDERS 3, 
Series A° 


Series B 


d 


PRISMS^ 

_ c 

Series A 


23.32 iT 0 * 04 ' 1 


21 { 1 + — — — > 
yl .08 


24.78 V 


■—0.021 23 {1+ 0.083 j. 


V 


0.5 


29.56 Y 


•-0.03 


2 ? { 1 + — — } 


V 


0.717 


a Cylinders of slenderness ratio 2 

^ Prisms of slenderness ratio 3 
c 

Concrete with river pebbles of average compressive 
strength 300 kg/cm 2 

d Concrete with granulated aggregates of average 
compressive strength 360 kg/cm2 


* ■ 
O' 

z 

UJ 

cr 

t- 

10 

UJ 

_J 

t/) 

z 

UJ 

H 

h- 

O 

LU 

a: 

Q 


M 1 -’' 





VOLUME V 


F iG-5-! DIRECT TENSILE STRENGTH OF CONCRETE IN CYLINDERS: DIA 
D = 5, 8 , 10. 12-5 AND I5CM OF CONSTANT SLENDERNESS H /D 2 ) 
RELATIVE TO THEIR VOLUME V (43; 



O I 2 3 4 S 6 7 89 SO 

VOLUME V 


FiG>2 DIRECT TENSILE STRENGTH OF CONCRETE IN 'PRISMS OF 
(SIDE LENGTH = 4 ,7, 10 » 12*5 AND 15 CM OF CONSTANT BLEND*. C 
NESS H/o-3) RE I ATI VE TO, THEIR VOLUME V 


(V 
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representative of the test data. However the form in 
which they specified viz, 

c = c 1 .. (5.3) 

where ?„nd k are constants 

implies that as the specimen volume increases the 
tensile strength decreases indefinitely. This is 
contrary to the fact that the strength should approach 
a limiting value as the size increases as has been 
discussed in the earlier lection 5.2. Consequently an 
attempt is made to fit the data to an alternative form 
of size mean strength relation given by 


= 0 


1 . 


{1 + 


K 


V 1 


n 


which has been studied in detail in the earlier chapter 
for the associated strength distribution function. 

The material parameters are evaluated by trial and error 
as follows. The criterion for the choice of the set 
of values of K and n is the closeness with which 

the theoretical strength distribution functions $ (<*) 
and the one obtained from experimental data agree with 
each other at various volumes. However no quantitative 
measures of closeness have been attempted and the 

rv r»nr* Af In r*P i & a 5=5 f* nl 1 nwfl . 



104 


Prom the Pigs. 5.1 and 5.2 depicting the experi- 
mental data relating the specimen size and the failure 
strength, the experimental strength distribution function 
$ (a) corresponding to various values of a is 
obtained from 

0) = — - .. (5.4) 

H+1 


where a = is the stress level 

n =number-cf specimens failed upto a stress 
level a 

IT = total number of test specimens. 

The theoretical distribution function $ (° ) is 

evaluated by using the result of chapter 4, as given 
by the Eqn. 4.22, viz, 


** (° ) - 


1- exp 


* 0 



for a > cr^ 


* 

for a < 


.. (5,5) 

The quantities n, a and K have been arrived at 
by trial and error so that the quantity $ {a) obtained 
from Eqns. 5.4 and 5.5 are in close agreement. The 
size-mean strength relation in terns of the material 
parameters a K and -n obtained by trial and error 
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as described above, in the case of cylinders (Series A 
end B), Prisms (Series A) are given in the last column 
of Table 5.1. 

5. THEORETICAL ATO OBSERVED DISTRIBUTION PUl'TCTIONS 

Pigs. 5.5 to 5.12 give the plots of $*(o) obtained 
c:ce orimontallv by using Eqn. 5.4- and that obtained using 
Eon. 5.5 with size-mean strength relation as given in 
bhe last column of Table 5.1. The Pigs. 5.5 to 5.12 also 
contain the theoretically computed failure frequency 
f (a) given by Eqn. 4.25 viz, 

y ( a ) = I { Ilstll > / ( ff _c i ) 1 /n-1 x 

emn [ -Y { 1 /n ^ _ a _ r ) Vn j 

iv. O’ -U 

for a > o T 

L .. (5.6) 

= 0 for a < 

Pig. 5.15 contains the failure frequency curves theoreti- 
cally evaluated from Eqn. 5.6, corresponding to various 
volumes in the cases of cylinders (Series A and B), The 
plots of $ ( 0 ) and f( 0 ) as given in Pigs. 5.5 to 5.15 
corresponding to various volumes of specimens are typical of 
the analytical study and as such only some of the results 
are presented. The plots indicate the increase in disper- 
sion in strength with decrease in specimen volume. 
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5.5 DISCUS SI Off 

m 

It can "be seen from the plots of $ (cr ) 

(Digs. 5.3 to 5.12) ; the strength distribution functions 
obtained theoretically and experimentally are in agree- 
ment with varying degrees at various specimen volumes. 
Taking for example the results of Series A cylinders, 
Fig. 5.3, qualitatively the plots of ( a ) are in 
better agreement corresponding to cylinder volume 

3 

V = 1.57 dm as compared to those in Fig. 5.4 and 5.5, 
indicating that the same values of the parameters 
a jj , F and n at various specimen sizes are inadequate 
to characterize the scatter with same degree of accuracy. 

It is interesting to compare the procedure 
adopted presently with that of Weibull as discussed in 
C hap ter 2, Section 2.25. In the Weibull’ s method of 
plotting the curves $ (a), the lowest strength a u 

is arrived at by trial and error so that the logarithmic 

A 

plot between log — ■ — and (a - a ) is nearly a 

1- $ (a ) u 

straight line. The values of the other parameters 
are fixed by the linear plot as discussed earlier in 
Section 2.2.3. Consequently in Weibull' s plotting 
only one parameter a u is fixed apriori by trial and 
error and the others are fixed by Weibull’ s plot. 

In contrast, in the reparometrized solution ns 
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need in the present charter, all the three parameters 
are obtained by trial and error. Even though while 
using a computer the difference in eomputational 
effort between the trial and err or procedure adopted 
presently and that in v fe lbull's graphical method may 
not he significant, Weibull’s graphical method is 
systematic in its application. However in the present 
case an attempt is made to utilize the results of the 
theoretical results of Chapters 3 and 4 and as such the 
results are applied to the direct tensile strength of 
concrete. However, if a more general strength distri- 
bution function, other than that of ■'Yeibull is sought, the 
associated parameters are to he found only by trial 
and error, unless some systematic procedure is developed 
depending on the nature of distribution in any specific 
application. 



CHAPTER SIX 


SIZJE EFFECT 0R_ TH E ER ROR AND RELIABILITY III 
PREDICT ! OH OF STREHGTH IN MAT ERIALS TESTING 


6.1 IHTRODU CTI Oil 

The importance of materials testing in the develop- 
ment of physical concepts and in engineering design can 
hardly he overemphasized. A true understanding and inter- 
pretation of the results of tests on materials is 
essential in assessing their behaviour and in design 

\ 

applications. Usually testing is so planned and 
motivated, that the results are representative of the 
material characteristics in a more general situation, 
while the tests should he simple to perform. Conseque- 
ntly, material testing is basically predictive in nature. 
In majority of the cases, there are several factors 
that are to he accounted for before accurate inference 
could he drawn from test results. The factors 
arise because of nonfeasibility of easily simulating 
all conditions of the prototype structure in the 
tests such as, the load, environmental and geometrical 
conditions that are present in a prototype or the 
original structure in which the material behaviour is 
to be predicted. The various factors affecting strength 
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rupture tests or Toy splitting cubes or cylinders. Each 
of these show different strengths. The kinetic effects 
on material testing are manifested by the variation 
in strength with rate of loading. As has already been 
discussed in Section 1.3 of the thesis the size effects 
on strength are essentially due to occurance of flaws 
and imperfections in materials. The flaw occurance being 
random, the strength can be specified only by a distri- 
bution function. In such' a case, it is of interest 
to assess the minimum number of test specimens that 
are to be tested so that the error in prediction of 
mean strength of the material does not exceed specified 
limits under different states of stress, for example 
flaws in concrete affect the tensile strength more than 
the compressive strength. Webster (120), studied the 
optimum sample size that corresponds to the minimum 
combination of testing as well as replacement cost, 
using information theory and Bayesian statistics. 

In the study of Webster, the specimen size is not a 
parameter that influences the cost. Since the number 
of specimens to be tested for predicting the mean strength 
is dependent on the specimen size adopted, an assessment 
of the optimum specimen size that corresponds to minimum 
cost of testing programme is of interest. The 
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aforementioned problems are discussed in detail in the 
following sections. 

6.2 SIZE AS A FACTOR IF SPECIFICATION OP MINIMUM NUMBER 
OP TEST SPECIMENS FOR GIVEN ERROR AND RELIABILITY 
IF PREDICTION OF STRENGTH 

The size effect on strength of various material and 
the origins of the same are discussed in Section 1.3. 

The general feature of size effect is that the mean 
strength and associated scatter decrease*' with increase 
in specimen size. Consequently a problem of interest 
is to find the minimum number of test samples that are 
to be tested §o that the population average strength 
can be predicted with a desired level of reliability 
and accuracy. Since the scatter or the coefficient of 
variation is affected by specimen size, the minimum 
number of test specimens to be tested also depends on 
specimen size. It is proposed to examine this problem 

V V 

making use of the test results of Kadlecek and Spetla 
( 43 ), on the tensile strength of concrete. 

Since the basic problem of interest is to correlate 
the sample and population averages, Student's f t' 
distribution may be used as was discussed by Cordon 
(ill), Beniamin and Cornell (112) and Hamilton (113) 
to correlate the error of the sample average, coefficient 
of variation and the number of tests as given by the 
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equation 

4- n O 

H * C-^ 1 ) .. (6.1) 


where 


and 


IT - number of tests 

t = value of Students ' t’ for N-1 degrees of 

freedom at a specified level of probability 
C v = coefficient of variation 
E = maximum percentage error of the sample 
ave rage 


The use of Student’s t distribution as given by Eqn. 6*1 
implies that the distribution of strength of the material 
is Gaussian with the coefficient of variation not known* 
Also it is assumed that the best estimate of the coeffi- 
cient of variation is the same as that given by the 
coefficient of variation of the sample. Results of 

V v 

Kadlecek and Spetla (43) from direct tension tests on 
concrete cylinders and prisms are shown in Table 6.1. 

The table gives the details of specimen sizes, average 
strength, coefficient of variation and the number IT 
of test specimens calculated according to Eqn. 6.1 
corresponding to an error E of 5% with 907. probability. 
’7h.il e calculating by Eqn. 6.1 for IT, the value of * t * 
is taken from the standard table of Student's distri- 
bution (111,113), so that the value of IT and t are 
eonsistant; i.e. the number of degrees of freedom R-1, 
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TABLE 6.1: HUMBER OB TEST SPECIMENS (H) REQUIRED TO BE TESTED 
TO ENSURE AH ERROR OB 5% WITH 90* PROBABILITY 
IH ESTIMATING THE MEM STRENGTH, AS MBECTED BY 
SIZE OB THE TEST SPECIMENS (DIRECT TENSION TESTS) 


Series A Series B 


* Average Coef.fi- Average Coeffi- 

,oxrme g-j_ ze cm stren- cient «■ stren- cient 
dm" 3 gth. p of vari- gth p of.vari- 

Eg/cnr ation% Eg/cm ation % 


CILINDERS 


0.24 

5 x 12 

22.8 

15.2 

26 

23.5 

12.6 

21 

0.81 

8 x 16 

23.4 

9.1 

1 1 

24.8 

10.6 

14 

1.57 

10 x 20 

23.2 

8.2 

9 

24.7 

7.8 

9 

3.07 

12. 5x 25 

22.2 

6.0 

6 

24.2 

8.1 

9 

5.30 

15 x 30 

21.5 

7.2 

7 

23.9 

6.7 

7 

10.60 

19 x 38 

21.3 

7.0 

7 

23.6 

5.5 

5 


P RISERS 


0.26 

** 

4x4x16 

30.0 

13.2 

21 

33.0 

10.4 

14 

1.03 

7x7x2 1 

29.6 

7.5 

8 

31.0 

9.1 

11 

3.00 

10x10x30 

28.4 

6.4 

7 

30.7 

6.7 

7 

5.85 

12.5x12.5 

x 37.5 

28 • 1 

5.7 

5 

28.9 

4.4 

4 

10.10 

15x15x45 

27.5 

6.1 

6 

29.6 

6.8 

7 


* diameter x length in the case of cylinders and side x side 
x length in the case of prisms 


** aspect ratio is different from the rest in the set 

Series A: Average compressive strength 300 Kg/cm 
Series B: Average compressive strength 360 Egcm 
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is so chosen that the value of h and the value of *t‘ 
from the Student's distribution table give a value of 
the coefficient of variation from Eqn. 6.1, very 
close to that given in fable 6.1 for various cases. 

It can be noted from the values of H calculated and 
given In Table 6.1, that the number of specimens to be 
tested for a given accuracy in prediction of strength 
varies with size. A smaller number of tests are necess- 
ary if larger specimens are used, than those required 
for smaller specimen sizes. Also the number varied 
with specimen shape and aggregate type, as a consequence 
of th^ change In the coefficient of variation. It is 
remarked that for a given specimen size, different number 
or sample tests are required in tension and compression 
Tor the saic error and reliability in prediction of 
mean strength for concrete because the flaw sensitivity 
of the material and consequently the coefficient of 
variation varies with state of stress. 

lor an illustration of the above remark the mini- 
mum number of specimens N, required to be tested to 
predict the mean compressive strength of concrete, 
calculated using the results of Rajendran (41), may 
bo used. Table 6.2 gives the details of test data of 
Rajendran and the values of N obtained by the application 
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of Student’s t distribution through. Eqn. 6.1, as "before. 


TABLE 6.2: EMBER OE TEST SPECIMENS (N) REQUIRED TO BE 
TESTED TO ENSURE AN ERROR OE 5% WITH 90% 
PROBABILITY IN ESTIMATING THE MEAN STRENGTH 
AS AEEECTED BY SIZE OE THE TEST SP ECIMENS^DIREC-T 
COMPRESSION TESTS) 


Volume 
Ou ” 

Size 

in. 

Average 

strength 

lb/sq* 

Coefficient 

of 

variation % 

N 

64 

4x4 

x 4 

4654 

2.936 

3 

216 

6x6 

x 6 

4450 

2.813 

3 

512 

8x8 

x 8 

4288 

2.799 

3 

100C 

10 xIO 

x 10 

4088 

2.714 

3 


It can he seen from Table 6.2 that the minimum 
number of test specimens did not differ corresponding 
to various specimen volumes as the decrease in coeffi- 
cient of variation with increase in specimen size is 
very insignificant. Also comparing the values of ’N’ 
from Table 6.1 and 6.2, it may be seen that more speci- 
mens are required to be tested in tension than in 
compression for the same error and reliability in predic- 
ting the corresponding mean strength. It should be 
noted, that this observation is by no means perfectly 
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conclusive since the material as well as the specimen 
volumes considered in Tables 6.1 and 6.2 corresponding 
to tnocic of Ta'llecek and Spetla (4-3) and Rajendran (41), 
are entirely different. 

Studies of Johnson (40) indicate that for a given 
mix of concrete, and for a given specimen size, the 
coefficient of variation is more in indirect tension 
tests than in compression tests. Consequently it may he 
inferred that more specimens are required to he tested 
n ‘n tension than in compression, to specify mean strength 
with a given reliability and error in prediction. 

however, since the state of stress in indirect 
tests is not purely tensile, the equalitjr of specimen 
volumes in compression and in ‘ indirect tests does not 
ensure perfect identity of test conditions. Since the 
material parameters in the strength distribution of a 
material are characteristic of the material, the coeffic- 
ient of variation under different stress states should be 
expressible in terns of the parameters. 

6., 3 OPTIMUM SPECIMEN SIZE AND SAMPLE NUMBER IN MATERIALS 

TESTING 

Since the mean strength and the scatter or coeffi- 
cient of variation of strength generally decrease with 
increase in specimen size, the number of test specimens 
N, to be tested to predict the mean strength of the 
material with a given error and reliability level varies 
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with the size of test specimens, As discussed in the 
preceding sections of the chapter, as the coefficient 
of variation decreases, the number of specimens llo he 
tested decreases; qualitatively, smaller the specimen 
size more the number to be tested. Apparently, in any 
programme of testing, one is confronted with a decision 
on the specimen size and the number of specimens so that 
the decision is optimal. The criterion of optimality 
could be v T ith respect to the total cost of testing 
which is the sum of the material cost and cost incurred 
in testing the specimens. Since, by adopting larger 
specimens (i.e, more material and more material cost), 
smaller number is required (i.e. less number of speci- 
mens and less testing cost), it is of interest to find 
the optimum sample number and specimen size so that the 
total cost of testing is minimized. This problem is 
examined in this section. 

The relation between minimum number of test speci- 
mens N and the coefficient of variation is given by 
Eqn. 6.1 (as discussed in Section 6.2), The coefficient 
of variation, in general decreases with specimen volume 
Y t The total cost of testing may be written as 

C T = ( K* V + 4 ) » 


. ( 6 . 2 ) 
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where = total cost of testing in Rs. 

K-j = cost of unit volume of the material in Rs. 

Y = specimen volume 

IT = number of specimens 

.* 

ana Kg = cost of testing a, single specimen. 

The quantity IT In Son. 6,2 is dependent of specimen volume 
Y and the coefficient of variation in strength C , 
through Eqn. 6.1. The possibility of an optimum speci- 
men volume Y existing is apparent, since by adopting 
specimens of large volume, the material cost increases, 
while the testing cost decreases; by adopting smaller 
specimens material ;ost decreases while the testing cost 
increases. The exact variation of total cost with 
specimen size V depends, however, on the values of K^ 
and Kg in Eqn. 6.2. 

Considering for example the direct tension testing 
of concrete, the number of test specimens IT, required 
to be tested to ensure an error of 5’/* , with a probability 
of 0.90 in estimating the mean strength is given in 

V ^ / \ 

Table 6.1, using the results of Kadlec eh and Spetla (43). 
To examine the total cost of testing as dependent on 
specimen volume adopted, the values of K 1 and Kg in 
Eqn, 6.2 are assumed and are as follows (for purposes 
of illustration). 
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= Rs. 0.5 per ctor of concrete , . (6.3) 

* 

-^2 = Rs* 5 per test specimen .. (6.4) 

so that 

= ( 0.5 V + 50 IT .. (6.5) 

In specifying the cost of material through K*, approxi- 
mate cost of T glue* necessary to fix the specimens to 
the end plates, so that the specimen could he tested in 
direct tension, is also included. Since the values of 
’IT’ for various specimen volumes ' V are alreadjr cal- 
culated according to Eqn. 6.1, and are given in Table 
6.1, the total cost may be evaluated using the values 
of and Kg as given in Eqns. 6.3 and 6.4. The 
total cost Crp according to Eqn. 6.5 is evaluated for 
the data in Table 6.1 (corresponding to series A and B 
cylinders and series A and B prisms)and the results 
are plotted in Bigs. 6.1 to 6.4. The plots show the 
variation of total cost in Rs, with specimen volume 
V(dm ). In plotting the Bigs. 6.1 and 6.2 the cost 

3 

corresponding to cylinder volume 0.24 dm Is not con- 
sidered: .-since' 'the’' cylinders have a different slenderness 
ratio from the rest o~f the specimens as is indicated 
by a double asterisk in Table 6.1. Similarly, the 

3 

cost of testing prisms of volume 0.26 dm is not 
considered in plotting Bigs. 6.3 and 6.4. The points 
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corresponding to actual calculated values of C^, for 
various volumes in the figures are joined by straight 
lines eventhough they could possibly lie on convex 
curves, as shown in dotted(profile only) if a consistant 
decrease in the coefficient of variation of strength is 
present with increase in specimen volume (as is theore- 
tically expected). The presence of local maxima, in 
the plots of the variation of with V in Pigs. 6*2 and 

6.3 is a consequence of the higher value of the coeffici- 
ent of variation of strength corresponding to that volume 
in the experimental results (possibly the limitations of 
accuracy of the experimental results). However, in case 
a consistant decrease in coefficient of variation of 
strength with increase in specimen volume is present 

(as is theoretically expected) such local maxima will 

not be present. Prom Pigs. 6.1 and 6.2 it can be seen 

3 

that the optimum volume of cylinders is 3.07 dm and the 

3 

optimum prism volume is 3.85 dim for the problem consid- 
ered using the results of Kadlecek and Septla. The 
corresponding sample sizes f N f are as given in Table 
6.1 corresponding to the optimum specimen volumes. 

6.4 DISCUSSION 

In this chapter some possible practical applications 
of brittle fracture statistics are discussed. While 
correlating the strength from tests to that in actual 
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structures, it is to “be noted that, especially in materi- 
als like concrete, quality control exercised in manufacture 
and placing can significantly affect strength as well as 
coefficient of variation (114 to 119). For example, the 
test results of Eadlecek and Spetla (43) and Rajendran 
(41) considered in the present study are for concrete 
manufactured and placed under laboratory conditions and 
consequently the quality control exercised can be expected 
to he good. However, the same concrete under field 
conditions may not reproduce the same scatter character- 
istics. nevertheless, the size effect on strength is 
prevalent in any case and the study carried out, to find 
an optimum specimen size for minimum cost of testing 
will be useful in planning a test programme economically. 

The basic principles underlying the example considered 
in the present study can be applied to other materials 
as well in standardizing the test programmes. Even though 
the assumption of normal distribution of strength in 
assessing the number of specimens applying Student's 
distribution is approximate (especially as applied to 
materials that exibit a significant skew in their 
strength distribution), it provides a simple means of 
fixing the number of test specimens economically in 
planning a test programme. 
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Through, the study of Johnson (40) indicating the 
existence of larger coefficient of variation in indirect 
tension testing of concrete than in compression (for 
same specimen size and mix quality) and from the present 
study as is discussed in this chapter based on the 
calculations of Tables 6.1 and 6,2, more specimens are 
required to be tested in tension than in compression to 
specify the mean strength with a given error and reliabi- 
lity in prediction. However } as already remarked, due 
to the presence of secondary stresses in indirect 
tension, the equality of specimen volumes does not 
ensure identity of test conditions in compression and 
indirect tension, further studies towards analytically 
obtaining the coefficients of variation under different 
stress states (like splitting of cylinders etc.) in terms 
of the basic material parameters, are necessary. Such 
studies in conjunction with those on the actual mechanisms 
of fracture involved under different states of stress in 
materials would give further insight in the statistical 
aspects of strength of materials as well as in consequent 
practical applications. 



CHAPTER SEVEN 


A PHMCMMCEDOG-ICAIi theory op fracture 
BEHA7I ' OT T OF, CONCRETE - LIKE MATERIALS 

7.1 INTRODUCTION 

In the preceding chapters of the thesis, the 
phenomenon of size effect on strength and associated 
scatter, is studied in some detail. Of related 
interest is the phenomenon of size effect on 
mechanical behaviour in materials. While size 
effects on fatigue strength and ductile-brittle 
transition in materials are known to be prevalent 
(58,15) as discussed in Section 1,4 of the thesis, 
systematic enquiries into the various aspects of 
size effects on mechanical behaviour of materials 
are scanty. This observation is somewhat surprising 
in the wake of the fact that size effect on strength 
has been so extensively studied and one could 
always conjecture logically that the same features 
of material structure like flaws, which affect 
strength could contribute in some measure to affect 
the other aspects of mechanical behaviour of 
materials as well. Dike size effect on strength, 
size effect on ductile-brittle transition, size 
effect on stiffness etc., are not encompassed by 
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classical theories of mechanical behaviour of 
materials. These phenomena cannot be rationally 
interpreted without resort to the study of inherent 
imperfections in materials. An insight into various 
aspects of size effects on stiffness, nonlinearity 
in st res s-st rain behaviour, ductile-brittle transi- 
tion etc., is essential in correlating the model 
and prototype tests, in addition to being of funda- 
mental value in understanding the behaviour of materi- 
als. Consequently an attempt is made in this thesis 
to study some aspects of this problem. 

In the following an attempt is made to propose 
3. unified theory that correlates the phenomena of 
size effects on strength and size effects on stiff- 
ness and ductility in materials. The theoretical 
development is applicable to materials in which a 
progressive breakdown in the internal structure of 
the material due to growth of microcracks is present. 
Also in the theoretical development, time, or temper- 
ature effects and multi axial stress conditions are 
not considered. In the formulation of the theory, 
some results of the theory of composite materials 
are invoked and consequently in the following sections 
a brief review of the background and implications 
of the results invoked is carried out for an appraisal 
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of the present contribution, 

7.2 THEORY OF C CMP OSI TE MATERIALS 

Though composite material systems like mortars 
and concrete are known to be in wide use from very 
early stages of technical activity, the increased 
use and innovation in several new composite materials 
like alloys, fibre reinforced plastics etc. of late 
has triggered active research into this area (101). 
"Thile the experimental and theoretical studies on 
the mechanical behaviour of composite materials are 
of interest in assessing their properties for 
engineering application, the subject is also of 
fundamental value since the related concepts provide 
a means of assessing the overall mechanical proper- 
ties of a material system in terms of its micro- 
structure. This aspect can be appreciated by observ- 
ing the gross behaviour of a polycrystalline aggre- 
gate. A polycrystalline material may be visualized 
as a complex material consisting of several randomly 
oriented grains, and knowing the properties of 
individual grains the behaviour of the aggregate may 
be assessed (100). The basic problem of interest 
in the theory of composite materials is to arrive 
at an estimate of the properties of the material 
in terms of those of the constituent materials and 
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the details of the geometry of the constituent 
phases. The following section briefly outlines 
methods of estimating the elastic moduli of a composite 
material, 

7 . 2 . 1 Methods oC Ev a luating the ELastic Moduli 

To estimate the elastic moduli of a composite 
material, the straight forward approach is to measure 
geometrically the relative locations and extents of 
the various constituent phases with known material 
properties and solve the field equations of elasticity 
in each phase and match solutions on the phase bound- 
aries. T ‘ r hile such an approach will be useful in 
studying materials where there is a definite pattern 
in phase disposition like in layered media, the same 
is practically impossible in the case of complex 
materials like concrete. In addition, in such materials 
the phase disposition is essentially random and for 
a formal approach to the study of such media, an 
appeal must be made to statistical continuum theories 
(102). In view of the fact that study of composite 
materials through a statistical theory is formidable 
in many cases, one is forced to go in for what are 
known as effective moduli of the composite material 
which are representative of the average behaviour 
of the composite material. Formally (100) the concept 
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of effective moduli can be given by 

"ij = X * J + 2 ”* ? lj ••• (7.1) 

where cr.^ = space average stress 

"ii = s P aCe average strain 
1 * 

X* and w = effective Lame’s parameters 
and A = space average dilatation 

Paul (103) and Hashin (104) derived bounds on the 
oCfective elastic moduli E and G using energy 
principles of elasticity theory, yhile fixing the 
bounds m elastic moduli by the energy principles, 
one is confronted with the task of finding an admissi- 
ble stress and strain field in the medium. To find 
the admissible stress and strain fields often the 
geometry of the phases in the composite materials is 
approximated by some simple geometrical model like 
snhere. Budianskg (105) attempted to evaluate the 
elastic moduli of a composite material with phase 
geometries being approximate spheres. The author used 
energy principles and self consistent approximation 
(106) to consider the interaction between the various 
phases. The results of Budiansty’s work for an ’IP 
phase medium are given by 



14T 


35T 

I 

i=1 


v i 


1 + P * - 1 ) 
G 


.. (7.2) 


I 

i=1 


v i 


1 + a K. 

H 

K 


- D 


.. (7.3) 


* 


a 


1 + v 
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2 (4-5 v ) ^ 

- T5 (TT/j - {7 - 4) 


where v^ = volume fraction of the i th phase 


G^ = shear modulus of the i th phase 

v = effective poisson's ratio of the 

composite material 

G* = effective shear modulus of the compo- 
site material 


K*»K, = effective Bulk modulus of the composite 

material* Bulk modulus of the ith phase 
' # respectively 

a ,13 = constants defined by Eq, 7,4 


IT = total number of phases 


In the limiting case of holes or voids in a two phase 
imtedium the effective moduli, K and G are given by 


( 105). 

0* = IU. zJL s l G (7.5) 

3 - v 

< 

K * = A ( . 1-2 Jl te l G .. (7.6) 

;v( 3-v) 



where v = volume fraction of voids 

G- = shear modulus of the solid phase 

Before concluding the discussion on the effective 
moduli of composite media, it may be mentioned that at 
present no sufficiently formal and simple approach is 
available to study the bulk behaviour of composite media. 
This may he noticed by the fact that the various bounds 
on affective moduli obtained by energy principles are 
too wide aoart to be useful (101,104), if one of the 
d has os happens to be voids or rigid. Of the various 
results available, Budiansby's (105) results are simple, 
event hough approximate and consequently the results 

given by Bans. 7.5 and 7.6 are used in subsequent 
development. 

7.3 A THEORY 0? FRACTURE BEHAVIOUR OF MATERIALS 

In the following, it is proposed to develop a 
unified theory that relates the size effects on strength 
with other aspects of mechanical behaviour like ductile- 
brittle transition, size effect on stiffness and non- 
linearity in stress-strain behaviour. The basic 
concept of the theory is that with increased loading 
microcracks and flaws develop and progressively decrease 
the stiffnc ss of the material as a consequence of 
progressive breakdown in internal structure. The extent 



of microcracking with loading is related semi-empiri- 
cally to the probability of failure at any stress 
level. The effective Young's modulus of the material 
in. the presence of microcracks is obtained by us in g 
Budiansky's (105) result given by Eqns. 7.5 and 7.6' 
for the effective elastic moduli of a two phase 
composite material in the presence of voids. The 
effective Young's modulus, so derived is shown to 
be dependent on stress level and probability of 
failure and is shovii to predict the nonlinearity in 
stress-strain behaviour, ductile-brittle transition, 
and stiffness as affected by specimen siae. The 
detailed development of the theory is shorn in Section" 
7 . 3 . 2 . , . 

7. 3. 1 Nature of Materials for which the theory is 
EBtfSa e 

The class of materials for which the theory is 
applicable is that in which permanent deformation due 
to microcracking takes place during sub critical crack 
growth. Concrete, rocks and fibre reinforced mortar 
are good examples of such materials, in which 
microcracks could be present in the material even 
before application of loading due to shrinkage and 
other effects. The microcracks progressively increase 
under loading and the microcracking- forms an energy 



dissipating mechanism other than surface energy. Such 
aforementioned materials which may he termed as 

,5 "I jfcs. mat eri al s 1 display a nonlinear stress- 
es rain "behaviour and behave in a semi-ductile manner 
in contrast to brittle materials in which the onset 
of crack growth and failure are synonymous. In the 
following theory the feature of subcritical crack 
growth due to microcracking is incorporated and hence 
the theory is applicable only to concrete like material 
in contrast to ideally brittle or ideally ductile 
materials. 

7.3.2 Development of the Theory 

As has been already indicated, the basic concept 
of the theory is to treat concrete like materials as 
composite materials with one of the phases being voids 
due to microcracking. It is essential that the 
microcracks be given due recognition since they signi- 
ficantly influence the overalls macroscopic behaviour 
of the material. The mioroc racks could be present 
even under no load due to differential shrinkage and 
consequent bond failures at the aggregate-mortar 
interface, or could develop during loading. In any 
case the microcracks propagate further with increase 
in applied load and influence the overall behaviour 
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and strength. The size effect on strength i,e, , the 
decrease in strength and associated scatter exist 
in concrete as has been already discussed in the 
preceding chapters. To keep the generality of the 
study, it is assumed that the strength distribution 
function characterizing the probability of failure 

•Wr 

( c ) , for a stress level less than or equal to 
c is given by 

* (e ) = 1 - exp {-v 2 ( CT )] .. (7.7) 

whore $ (cr ) = probability of failure when the 

stress lies any where between 0 and a 

V = specimen size (volume) 

f (a) = is any function of a that increases 

with o and is such that it satisfies 
suitable conditions to make , (a), 

o, distribution function (? (a) could 
he, for example ( cr - cfa ) m as in 

,r . r e ibull f s di o t ributiofi ) . 

It could be observed from Eqn. 7.7 that the probabi- 
lity of failure <> (c) increases with increase in 

stress level for a given specimen size and the proba- 
bility of failure increases with specimen size for 
a given stress level. The specific form of ? (o ) 
is not necessary in"the following theory, 

A typical test on concrete in compression with 
measurements for linear, lateral and volumetric 
strains, fig. 7.1 reveals that the stiffness of the.. 
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material goes on progressively decreasing with applied 
load and the stress - strain curve is nonlinear at 
each stage of loading. The nonlinearity in the stress 
strain behaviour is attributable to the irreversible 
microcracking (80), the initiation of which becomes 
significantly ’manifest at applied loads of above 
30% of the ultimate 1' ad. With further loading, the 
stiffness progressively decreases and when the applied 
load Is about 70 to 90%, (107) of the ultimate, 
considerable breakdown in the internal structure 
takes place when the well developed microcracks 
start joining and form fissures so that at the ulti- 
mate stress level, the material system tends to 
form a mechanism, Further, an examination of the vari 
at 1 on of vaLumetri ? strain with load indicates a 
progressive decrease in magnitude because the system 
is in compression, but with the onset of disconti- 
nuity in the material i.e., at about 70% of the 
ultimate load, the volume change reverses in nature. 
While many attempts have been made earlier to express 
the stress-strain relation of concrete by various 
analytical expressions (108), none of these are based 
on any formal theory. 1 formal approach to the 
problem is however, extremely difficult, but it 
appears that the same may be approached in a rational 
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manner with some suitable approximations and simpli- 
fications as discussed below. 

In a study, wherein a theoretical evaluation of 
the stiffness of concrete as dependent on applied 
stress level is attempted, the following points should 
be incorporated. While a progressive decrease in 
stiffness of the material with increasing stress indi- 
cates that there is an increase in the volume of 
microcracks or voids in the material, the decrease in 
the gross volume change initially in the case of 
compression tests indicates, Fig. 7.1, that the solid 
phase decreases in volume with increased compressive 
stress. The macroscopic ally observed gr-'ss volume 
change is the net effect of increase in microcrack 
or void volume and that due to decrease in the solid 
phase. However at' stress levels beyond about 70% 
of the ultimate load, the nature of variation of 
volumetric strain with load changes, as shorn in Fig. 7.1, 
indicating that at the loads nearing ultimate load, 
the increase in volume of microcracks or voids dominates 

the decrease in volume of the solid phase. In the 
case of tension tests, the phenomenon can be expected 
to be quite anologous to what is presently discussed, 
except that the gross volume increases monotfaiically 
up to failure. 
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As was mentioned earlier, Eqn. 7.7, the stress 
dependent risk of failure of specimen of size V can 
he given hy a distribution function of the form, 

* (o) = 1 - exp { -Y f (a) } .. (7.8) 

Since concrete-like materials display a progressive 
breakdown in internal structure with increased stress 
a , it can be conceived that various elements consti- 
tuting the material progressively fail and contribute 
to the microcrack or void volume (109). Consequently 
tie materials do n n t obey the weakest link concept. 

!Po quantify the phenomenon of progressive breakdown 
of internal structure the cumulative probability of 
failure $ (a) can be visualized as follows. If the 

material of volume Y can be divided into IT elements 
and n( a ) elements are likely to fa£Lat an applied stress 
level o , then the cumulative probability of failure 
(<*) may be written as (109) 

(c i) = 1 - exo { -Y f (a ) > = ^ • (7.9) 

v ' I 

or n (a) [ 1 — exp { — Y f ( ® ) 3 3 «• (7.10) 

If n (c) elements get ruptured during loading upto 
a stress level a , then the volume v Q occupied by 
the microcracks can be given by 
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v Q = 1<" 0 n (cr ) = > 0 1'T [1 - exp {-V f (cr ) >3 . . (7.11) 

where k = a constant that denotes the mean volume 
of a single microcrack. 

The volume occupied by the n microcracks "being 

v , the effective volume of microcracks or voids that 
o’ 

contributes to the loss in stiffness can he given by 

v = k 1 v Q ..(7.12) 

where v = effective volume of microcracks or voids _ 

that decrease the stiffness of the material 

and k^ = is a constant 

It could be seen that Eqns. 7.9 through 7.12 correlate 

the probability of failure with effective void volume 

in the material. With this estimate of void volume, 

one may find the effective Young' s modulus E in the 

presence of microcracks by using Budian^.fey' s (105) 

results, given by Eqns. 7.2 and 7.3* Budiansky' s 

* * 

values of effective bulk and shear moduli K and G , 
respectively can be rewritten as 

= IIIjzJLI & .. (7.13) 

3 - o 


_ 4(1 - 2C) (1-fi_ 

c (3 - c ) 


.. (7.14) 


where c 

G 


effective void or microcrack volume 
fraction - 

V 

the shear modulus of the solid phase 
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The effective Young' s modulus E can he obtained from 
Eqns. 7.13 and 7.14 and is given by 


_ 9 K G _ 56(l-2c) (1-5) 

3K +G (3-c)c +4(l-c)(3-c) 


G .. (7.15) 


or 

E * = 36 ( 1-2c) (1-c) . _E 

( 3-c)c +4( 1-c)( 3~c ) 2( i+v ) 


.. (7.16) 


2i±lz2cLIlzcl (_JL_) 

(3-c)c +4 (i-c)(3-c) * 


.. (7.17) 


where v =0.2, the Poission's ratio of the solid 
phase 

and E = the Young' s modulus of the solid phase 

It is fortuitous that v = 0.2 used in Eqn. 7.17, 
following the analysis of Budiansky, happens to be 

a good approximation for the Poisson's ratio of 
concrete during the initial stages of loading. 

It can be seen from Eqn. 7.17 that as c^ -fc^e 
fractional volume of voids, increases E decreases and 
approaches zero when c equals a value 0.5, implying 
that E reduces to zero or failure of the material 
occurs, when the effective void volume reaches 50% 
of that of the total material. This is contrary to 
usual observations and. is an outcome of the theory 
of Budiansky (105), as remaked by he himself. 
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Concrete, for example has void or microcrack volume 
at failure much less than 50% of the total volume. 

This is essentially the reason why the constant k.j 
was introduced empirically in Eqn. 7.12. The constant 
k,j ur^viOos for the necessary magnification of the 
actual ,r oia volume so that Budiansky's theory is made 
suitable and may be evaluated by experimentally 
measuring the fractional void volume at failure 
and setting 

V hauus = 0.5 ..(7.18) 

T 

failure = void volume at failure 
V = total volume 

ft 

The effective Young's modulus E as given by Budiansky's 
theory, Eqn. 7.17 can be evaluated if c, the 
fractional void volume is known. However, the 
fractional void volume j dependent on stress level 

ft 

a , probability of failure $ (a) and specimen size 

Y, is already given by Eqns. 7.11 and 7.12 and is 
given by 

** "V* 

c = |*-12 S [1- exp {-V f(a)}] ..(7.19) 

ft 

The effective Young's modulus B is given by 

E * , 36(1-2c) (1-c) ( E } = 15 E(1-2c)(1-H) 

( 3-c)c + 4 ( 1-c)( 5-c) (5-c)c+4(1-c)(5-c) 

.. (7.20) 
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whore o is given by Eqn. 7.19. 

In the -iqn. 7.20 above, the effective Young’s modulus 
* 

E of the material is expressed in terms of c, the 
fractional void or microcrack volume. Since c is a 
function of applied stress level a and the specimen 
sloe 7, the Young's modulus E* is a function of specimen 
sine end stress level. Consequently the theory 
developed yields a size and stress dependent value of # 
Young’s modulus E for materials in which microcracking 
ir prevalent, further implications of Eqns. 7.19 and 
7.70 will be discussed in detail in the following 
c h M o ter. 

f.4 DISCUSSION 

In this chapter an expression for the effective 
Young’s modulus of a material in the presence of 
microcracking is developed using an estimate of the 
volume of microcracks as dependent on stress level 
and specimen size. The growth of microcracks is assumed 
to contribute to the effective void volume in the 
material. The material is considered as a two phase 
composite material in which one of the phases is voids 
and the related result of Budianslcy (105) for effective 

Young* s modulus is invoked in the theory. The choice 
of Budiansky's result is purely for its simplicity 
and the logical basis on which it is developed. 
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Hovrever, to compensate for the fact that the fractional 
void volume at failure is much less than that 
predicted by Budiansby's result? an additional constant 
lc.| was to be introduced in the theory which could 
bn avoided if a more general and accurate result for 
effective Young's modulus of a composite material 
in the presence of voids is available. The impli- 
cations of the theory developed in terms of the 
nonlinearity in stress-strain behaviour, size effect 
on stiffness and ductile- brittle transition will be 
discussed in detail in the following chapter. 



CHAP TER EIGHT 

THEORY 0? SIZE EF FECTS ON MECHANICA L 
BEKAVKUR QE CONCRETE-LIKE MATERIALS 

3.1 NONIINEAKETY IN STRESS'* SI RAIN BEHAVI OJR 

Truesdell (66) remarks that the various theories 

of nonlinear constitutive laws can he classified to 

fall into two groups depending on the basic approaches 

In either of the aporoaches, a description of the 

quality of the material and a method of quantifying 

it are attempted. In the classical theory of 

materials vis, elasticity, _ the quality of material 

(perfect elasticity) is quantified in the simplest 

manner by linearity. Starting from this basis, in the 

first approach, additional qualities to the material 

/ 

are assigned while the quantification is still linear 
as in plasticity theory; in the second approach, the 
quality, namely elasticity is preserved but the quanti 
fying measure is changed as in the theories of finite 
elasticity. The suitability of the two approaches, 
however is governed by the material for which it is 
applied; ’Viewing in this light, the following 
discussion of the results of the preceding chapter 
can be observed to fall into the first category in 
the broad classification of Truesdell. 
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8.1.1 Nonline arity due to ^Progressive Breakdown 
i n Internal Stmctirrp 

The Eqns. 7.19 and 7.20 relating the effective 
modulus E and effective fractional void volume c, 
derived in the last chapter are as follows: 

c = k 1 k Q N ** (a) = k., k Q N [1- exp {-Vf (a)}] 

.. ( 8 . 1 ) 



(3-c)c+4(l-^)(3-c) 


.. ( 8 , 2 ) 


The stress-strain law of the material as described 
by the above Eons. 8.1 and 8.2 is essentially non- 
linear, since the elastic modulus E contains a 
through c. It is obvious that the nonlinearity is 
an outcome of the fact that the microcracks which 
contribute to the void volume, increase with stress 
0 . Erom Eqn. 8.1, it can be seen that as a increases, 
c the effective fractional void volume increasea and 
as c increases E* decreases, as given by Eqn. 8.2. 

When c tends to a value 0.5, E reduces to zero. The 
variations of c and E with 0 and the stress -strain 
relationship according to Eqn. 8.2 are shown schemati- 
cally in Pigs. 8.1 to 8,3. The slope of the stress- 
strain curve given in Pig. 8.3 denotes E , the 
effective Young's modulus. The predictions are 
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*1 ( c ) = 1 " ®XP { - ^ f(a) } .. (8.3) 

$ 2 = 1 “ { - V 2 f (a ) > .. (8.4) 

*1 (o) ^ (o ) .. (8.5) 

^ 1 

* 

1 ( c ) = the probability of failure for any 

stress level less than or equal to 
o , of specimen of size 

* 

9 2 ( a ) = the probability of failure for any 

stress level less than or equal to 
a , of speciment of size Yg 


Ion si do ring that the number of macroscopic volume 
elements in the material per unit volume to be given 
by n , the total number of volume elements N in 
Son. 8.1 is given by 

IT 1 = n 
hg = n 

size Vg ..(8.7) 

whore 

N,, = total number of macroscopic volume elements 

in the specimen of size 1/t 

Up = total number of macroscopic volume elements 
d in the specimen of size Yg 

The void volumes at any stress level a in either of 

the two specimens of sizes Y-j and Yg can be obtained 

relating the probabilities of failure (o) and 

( 0 ) to the number of ruptured elements in the 

two specimens, as follows: 


Y.| in the case of specimen of size Y 
Vn in the case of specimen of 


. .( 8 . 6 ) 
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^ (a) 

^(a) 

= -TT .. (8.8) 

n 7, 

1 

and 

* (a) 

2 

n 2 (a) , 

(8.9) 

n Y 2 

where 

* 

n 

= number of macrocropic volume 
elements per unit volume 


n 1 ( a ) 

= number of elements that get ruptured 
at a stress level a in the specimen 
of size 7,j 


CvJ 

3 

= number of elements that get ruptured 
at a stress level a in the specimen 
of size Y 2 


Since the raptured elements are considered to contri- 
bute to the effective microcrack or void volume, 
as discussed in the previous chapter 


v 1 = k 1 k 0 n-jCa) 


- k i 

k o n 

V 1 

$ 

1 

(a) 




[1- 


= k 1 

k Q n 

V 1 

- exp 

a$i4 similarly 





* 


[ 1 


V 2 = 

k o 11 

- exp 


- Y 1 f (o)}] .. (8.10) 

{- Y 2 1 (a)}] .. (8.11) 


where v, - effective volume of microcracks or 

1 voids in specimen of size 7^ 

Vp = effective volume of microcracks or 

2 voids in specimen of size 7 2 

k * constant denoting the mean volume of 
0 a single microcrack (defined in the 
previous chapter) 
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a,na lc 1 - semiempirical constant defined by Eqn. 

7.18 in the previous chapter 

i! rom Eons. 8.10 and 8,11, the fractional effective 

void volume in either of the specimens can be written 


as 

? i 

, , * 



{-Y 1 f(a)}] .. (8.12) 

C 1 = 

7 T 

= k. k n 

f 0 

[1- 

exp 

and 








# 



{-V 2 f(a ) }] .. (8.13) 

O 

ro 

li 

1 o 

t*. 

■ k 1 k o n 

[1- 

exp 


Prom Eons , 8,12 and 8.13, for a given a , since 


2 


\ *r 


‘ 1 


or) / 


.. (8.14) 


The Eqn. 8.14 implies that, for a given stress level 
c , the effective fractional void volume is more 
for a larger specimen. Prom Eqns.8.2 and 8,14, it can 
be observed that the effective Young's modulus E , is 
less for the larger specimen, at a given stress level 
a . If E denotes the slope of the stress-strain 
curve, then the stress-strain curves corresponding to 
specimens of different sizes would be, according to 
Eons. 8.2 and 8.14, as shown schematically in Pig. 

8.6. The Pig. 8.6 shows the predicted size effect 
on stress-strain behaviour. The mean failure stress 
is denoted by E reaching a value zero. The larger 
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specimen has a lower failure stress as E* approaches 
zero more rapidly with a according to Eqns. 8.2 and 
8.13. 

8.8 SIZE EFFECT ON STIFFNESS 

Fig. 8.5, which shows schematically the size 
effect on stress-strain behaviour also shows the 
si •’ o effect on stiffness. From the discussion in the 
preceding section and the development of Eqns. 8.2 
to 3.14, it follows that, at any given stress level, 
the larger specimen has a. greater effective void 
volume and hence a lower effective Young’s modulus 
E , implying that larger specimens are less stiff 
as compared to smaller specimens. 

0.3 SIZE EFFECT ON DUCTILE-BRITTLE TRANSITION 

The effects of temperature and state of stress 
on ductility of materials are well known. However, 
the effect of specimen size on ductility is less known 
and is essentially brought to light by Glucklich (58, 
61,62). In the following a brief account of the 
observations and theory of Glucklich is given. The 
predictions of the present theory with respect to 
size effect on ductile-brittle transition are compared 
with those of GluclcLich in Section 8.3.2. 
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8 . 3 . 1 jgkxU er view Point of Strain Energy -Size Effect 
GlucKLich (62) examined the load deflection 
profiles of notched beams of different sizes and the 
stress-strain curves as obtained from notched beams 
gyp sum mortar loaded by a spring. In addition, 
GlucKLich studied the effect of loading with a spring 
in series in compression and tension on specimens of 
Gypsum mortar. Glucklich noticed that the load 
deflection profiles of notched beams depicted that 
greater energy was absorbed before failure in smaller 
specimens. The presence of a spring in series with 
a notched beam had the similar effect yiz, specimens 
with no spring displayed greater ductility and strength 
than those with a spring. The observations indicated 
that the presence of a spring in series with a 
specimen under load is essentially same as having a 
larger specimen with reference to strength and 
ductility. Larger specimens showed less strength 
and ductility as compared to smaller specimens. So 
also, comparing the behaviour of specimens with 
identical size but with and without a spring in series, 
the specimen tested with a spring in series showed 
less strength and lower ductility. The above phenomena 
are commonly termed as Strain Energy Size Effect, 
by Gluchlich (58,62), The theory advanced by 



Glucklich to interpret these phenomena is essentially 
based on the total energy of the system under loading 
and balance of energy during crack propagation. The 
theory of GluckLich, briefly, is as follows. In an 
ideal Griffith material, the surface energy is the 
n nly mechanism of energy dissipation and a single 
potential crack and its propagation is considered. Also 
the ideal brittle material is homogeneous with a uni- 
form energy demand during crack propagation. Conse- 
quently an equilibrium exists between crack length 
and external load during stable crack propagation. 

As the load and crack length increase, a stage is 
reached when the strain energy release rate exceeds tiiat 
of energy absorption due to surface energy and sponta- 
neous fracture occurs. However in materials like 
concrete and many other real materials, the energy 
balance is not as described above because of the 
presence of microcracking and the presence of harder 
phases. Whenever a propagating crack encounters a 
harder phase greater energy is required to force the 
crack further. Also when once the crack encounters a 
weaker zone, the excess energy, if present w^rks 
against the uncracked material. Consequently in such 
materials the total energy of the system which acts 
as a reservoir determines the mode of crack propagation 
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and failure. If a system has higher energy level due 
to one rresence of a spring in series or a larger 
specimen, the excess energy drives the crack through the 
harder phases and the following kinetic energy advances 
fracture, favouring the conditions for unstable crack 
propagation. Because the deformation is smaller during 
unstable crack propagation, the presence of the energy 
reservoir due to spring or a larger specimen size 
brings down the ductility of the system also. Based on 
the above concepts, C-luckLich proposed a transition 
size anologous to transition temperature, which is 
schematically shown in Big. 8.7. 

8.3.2 Evaluation of the Present Theory in the Light 
of the Concept of St fain Energy Size Effect 

Taking total area under the stress-strain curve 
upto to the failure stress, i.e. the energy absorbed 
before rupture as a measure of ductility it can be 
noted that the size dependent stress-strain curve shown 
in Eig. 8.6 and discussed in Section 8.1.2, 
explains that larger specimens display less ductility 
as compared to smaller specimens. This feature is 
Schematically shown in Fig. 8.8. It should however be 
noted that the theory is applicable upto maximum or 
failure stress, but not beyond, for example, in the 
strain softening range. The viewpoint advanced in 
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the present theory? namely, that the microc racks c-ontri— 
"bute to void growth which lowers the material stiffness, 
is different from that of G-luckLich wherein the total 
energy content of the system and the energy balance, 
determines the fracture mode. In the present theory, 
else effect on strength, stiffness and nonlinearity 
in material behaviour are inherent in the material 
that displays size effect on ductile-brittle transition. 
While the theory of GluckLich emphasises on the basic 
mechanism of fracture in the materials, the present 
theory examines the behaviour essentially on a macro- 
scopic level and relies on the quantification of 
fracture process on the average, incorporating various 
known characteristics of deformation of the material 
like microcracking and volume change within the frame- 
work of the theory of composite materials. Both 
theories being qualitative at present, a detailed 
comparison is not possible. A more specific comparison 
needs further study of both the theories as applied 
to specific materials, 

8,4 E3CPERIMMTAL EVIDENCE IN SUPPORT OF THE THEORY 
Systematic experimental investigations into the 
size effects on strength, stiffness and stress-strain 
curve together on a^material are not reported in the 
literature. A preliminary experimental investigation 



169 


carried out on fibre reinforced cement mortar specim- 
ens, to examine the predictions of the theory is 
described helow. fibre reinforced cement mortar 
is a semiductile material which has stable crack 
p r on an; at ion and admits a progressive break down in 
internal structure under increased loads due to 
microcracking in mortar as well as the bond rupture 
at fiber - mortar interface. Since the strains are 
relatively larger, the experimental measurements are 
more accurate as compared to concrete and similar 
materials in tension. 

8.4*1 Details oT Tests on Steel fiber Reinforced 
Cemen t Mortar "Br j qu ett e s 

Pour series Of tests on fibre-reinforced cement 
mortar briquettes in tension were carried out with 
briquette thickness as the variable. The cement 
mortar prepared was of 1:3 mix consisting of sand 
of fineness modulus 9.1 and standard P ortland C ement. 
The embedded steel fibres were of 18 gage and yield 
strength 38 Kg. In series 1,2 the fibre content was 
8% by volume while in the series 3 and.'! it was 16% 
by volume. In each of the series the briquettes 
were prepared and tested under identical conditions. 
The briquette sizes are shown in Pigs. 8.9 and 8*10, 
Table 8.1 gives the ddtails and number of specimens 
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TABLE 8.1: DETAILS 0? TEST SPECIMENS 0E STEEL FIBER 
REINFORCED CME1TT MORTAR 


Series 

Fibre 

volume 

% 

Briquette 

thickness 

inches 

Number of 
specimens 
tested 

1 

8 

0.5 

3 

2 

8 

0.75 

3 

3 

16 

0.5 

3 


4 


16 


0.75 


3 
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STRAIN CURVES OF FIBRE REINFORCED 
WITH SPECIMENS OF DIFFERED SlZtS 
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tested in each of the series. 

8.4.2 Test Results 

The stress-strain curves as obtained from tests 
with ori riiiet fee thickness as the variable are shown 
in Pigs. 8.9 and 8.10. The Pigs. 8.9 and 8.10 clearly 
indicate that a size effect on stress-strain behaviour 
and strength exists. Also smaller specimens have 
hi -her strength, larger ductility and greater stiff - 
ness. The strains plotted in the Pigs. 8.9 and 8.10 
are the average values of the three specimens tested 
in each case and consequently indicate the average 
behaviour. 

8.5 C.QMPARISOJT OP EXPERIMENTAL OBSERVATIONS AND 
THEORETIC AL PREDICTIONS 

The theory of size effects on materials proposed 
in the thesis predicted smaller stiffness and less 
ductility associated with larger specimens. The 
predictions are qualitative and are in agreement 
with the experimental observations on fibre-reinforced 
cement mortar briquettes in tension as shown in 
Pig. 8.9 and 8.10. The theory implied that the failure 
of material and the effective Young’s modulus reaching 
zero are synonymous. Further, the descending portion 
of the stress strain curve, i.e., the strain soften- 
ing portion is not included in the theory discussed 
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and the predictions are valid only upto the maximum 
stress level. However, the Figs. 8.9 and 8.10 indicate 
higher ductility of smaller specimens, even by con- 
sidering the area - under complete stress strain 
curve as a measure of ductility. 

8.6 DISCUSSION 

Ihe expression for effective Young 1 s modulus 
of material as given by Eqns. 8.1 and 8.2, predicted 
a '<'"U vtely the size effect on stiffness, stress-strain 
behaviour and ductility at failure as evidenced by the 
experimental results discussed in Section 8.4. Ihe 
predictions were essentially qualitative since the 
constants , k Q and N are to be fixed suitable to 
any particular material. Ihe constant k Q denotes 
the mean volume of a typical microcrack, k^ is a 
constant that was empirically introduced to transform 
the microcrack volume to the effective void volume 
at failure as was discussed in Section 7.3.2, Eqns. 

7,12 and 7.18> while K denotes the number of elements 
into which the material is devided in the development 
of the theory and is only of analytical significance. 
•These constants may be considered as parameters 
characterizing the property of the material for which 
the theory is applied. While at present no specific 
suggestion can be offered as to the experimental 



detexmination of the constants, the quantity ’ k.j k Q N* 
ir.-r/ to taken together to denote a single constant of 
the material and the same may he quantified in any 
specific application. The qualitative implications 
of the theory namely, the greater stiffness, greater 
ductility associated with smaller specimens are 
however valid for any positive practical values of 
the constant k 0 K. 

The theory is not very rigorous in its formu- 
lation and is developed hy fitting several physical 
'•'"'nr epts within a logical frame work. Concrete-like 
priori -iis are ‘oasidered as two phase composite 
materials with one of the phases as voids. The micro- 
Ci-'.chs in concrete like materials are related to the 
voids in the cornu osite material. It is not implied 
that the correspondence between microcracks and voids 
is by any means exact; nevertheless, the correspond- 
ence Is invoked only to quantify the stress depend- 
ence of Young's modulus in such materials. The proba- 
bility of failure which is dependent on stress level 
and specimen size, assumed to be known, is related to 
the microcrack volume and hence to the effective 
Young's modulus of the material. The theory so 
developed adequately explained the size effect on 
stiffness , stress-strain behaviour and ductile-brittle 
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transition in concrete like materials. The transition 
from ductile to brittle as the specimen size increases 
is essentially gradual. The predictions of the theory 
arc observed to be oualifcatively in agreement with 
observations in tests on fibre reinforced cement 
mortar briquettes, further exhaustive tests are 
essential, especially with measurements on volume 
changes and extent of micro cracking at various stress 
levels, for a quantitative evaluation of the theory. 



CHAPTER NINE 


I.IATERIA1 CHOICE EOR FRACTURE RESISTANT 
DESfoiTW pressure vessels 


9 . 1 INTRODUCTION 

In the previous chapters, statistical aspects of 
size effects on strength and fracture behaviour have 
been investigated. In this chapter a study is made on an 
aspect of fracture resistant design as an illustrative 
example of use of fracture mechanics in engineering design. 
The presence of crack like defects in structures like 
pressure vessels, pipe lines, ship hulls etc., can lead 
to catastrophic failures, if not g-urded against in the 
design process. The conventional design process based 
on specified nominal stress level is inadequate unless 
the absence of cracks in structures is guaranteed during 
fabrication or subsequent service life. However, cracks 
in structures can be initiated and propagated during 
the fabrication as well as duringtteir service life. 

In the presence of a crack the integrity of a structure 
is impaired depending on the severity of the crack 
with reference to its location, size ana shape; Appli- 
cation of fracture mechanics concepts in structural 
design enables the engineer to assess the tolerable cracks 
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in a structure working at a particular stress level or, 
conversely, the allowable stress or load level in the 
presence of a crack of known characteristics. In the 
actual design problem where there is a likelihood of 
brittle fracture oc curing, the fracture toughness and 
crack propagation characteristics of a material to be 
used are to be considered in the design process in 
addition to the other properties of the material like 
unit cost, density etc. It is attempted in this chapter 
to cjevelop a design procedure for fracture resistant 
design an applied to pressure vessels. The emphasis in 
the design procedure is on the optimal material choice 
from an available set of materials with known properties. 
T.n the foil owing, the need for fracture resistant design 
in pressure vessels is briefly discussed. Subsequently, 
a specific problem of pressure vessel design in the 
presence of a part-through crack, with optimal material 
choice is discussed in detail. 

9.2 DESIGN CRITERIA FOR PRESSURE YESSELS 

Several failures of welded pressure vessels, oil 
and gasoline storage tanks and pipe lines (121,, warrant 
the need for fracture resistant design of such structures. 
While it is known that crack like defects originated 
either due to defective welding or during service life 
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are Imown to initiate brittle failures, the exact 
origin oi the strength impairing defects is not clearly 
known. However, the presence of cracks in structures 
can initiate fracture £ ' j ven at service stresses much 
lower than the yield stress of the materials. The 
conditions for fracture in structures like pressure 
vessels and pipe lines are more favourable for brittle 
fracture oocurance where the inherent ductility of the 
material is impaired due to various reasons. In the 
case of pressure vessels operating in a nuclear atmos- 
phere, the neutron bombardment on the vessel material 
improves the yield and ultimate strengths of the material 
but significantly lowers the ductility and impact 
strength (122). Also, in the manufacture of pressure 
vessels, plates which are initially flat are bent into 
desired shape resulting in the cold forming of the 
material and as a consequence the material is prestrained 
which raises the yield point hut reduces the ductility. 

In such cases the properties of the material prior to 
deformation need not be representative of those in the 
actual structure end this is more severe for fracture 
resistance (122). Pressure vessels under external shock 
loads and internal pressure cycles must have adequate 
ductility to operate economically at service temperatures. 
In a design attempted to realize such an end, the material 
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choice for the pressure vessel must he made with adequate 
fracture resistance. 

9.3 PRACTURE RESISTANT DESIGN OP PRESSURE YESSEDS 
9*5.1 Problem of fracture Resistant Design 

Wliile designing structures applying the concepts 
of fracture mechanics, the existence of cracks in 
structures is inherently admitted. In the design 

process one assesses the allowable stress level and 
the associated service life in the presence of flaws. 

The most important factor in the design is an estimate 
of the initial crack size, location and shape. Unfort- 
unately the initial cracks often are very small, that 
they cannot be directly measured easily or could be 
observed by any nondestructive testing techniques. 
Consequently an estimate of the initial crack size 
and shape could only be mate by an intelligent guess based 
on experience, observation of the previous designs and 
some knowledge of the fabrication and welding procedures 
adopted. Even though no formal procedure exists for 
initial crack size estimation, proof testing can be 
of some help in many cases (121). The method of 
assessing initial crack size using proof testing 

will he discussed in the problem attempted subsequently 
In this chapter. In designing pressure vessels for 
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fracture resistance, two approaches can be followed 
(i) fail safe design (ii) Safe life design. The most 
critical crocks in pressure vessels are those which 
initially penetrate through the vessel wall only partly, 
called ” part -through cracks” . Those cracks which are 
tirron.gh the wall thickness fully are usually detected 
by a leakage and could atleast be noticed. In contrast, 
the part-through cracks even though of such size initially 
that they are harmless, could grow during service life and 
may finally be of such size that catastrophic fracture 
occurs. In the fail safe design, the material and the 
wall thickness are so chosen as to ensure that the crack 
cannot catastrophic ally nrop agate unless it grows to such 
size that it penetrates completely through the wall 
thickness so that a leakage is indicated before fracture. 
In the safe life design, for a chosen material and wall 
thickness, the service life is estimated before fracture 
occurs and the structure is put out of service before 
the end of the design or actual life. The service could 
he resinned if the grown crack is of such size that it 
could he located and accurately measured and if it is 
found that further life can he ensured before rupture. 

It is the aim of the present chapter to suggest & design 
method so that the material choice could he optimally 
made for the design of pressure vessels hy either of 



ttie approaches viz, the fail safe and safe life designs* 
The method is illustrated "by a specific example in the 
following sections. 

3»3»2 All_^a ilple_ of the Design of a P ressure Vessel 
j-J L-tiie. P re sene e of a part -through Crack 

A circular cylindrical pressure vessel of length 

1000 Inches with an internal pressure of 2000 psi. is to 

he designed, i.e. , the material and thickness of the 

shell, are to be chosen. Both safe life and fail safe 

designs are to be attempted and it is given that the 

part-through crack is semielliptical in shape and the 

ratio of flaw dent! to the major axis is 0.4 as shown 

in fig. 9.1. The nominal working stress should not 

exceed sixty percent of the yield stress of the material 

used. It is required to find an optimal material 

(i) When the vessel diameter is fixed at 50 u (it turns 
out to he safe life design as will be shown below) 

(ii) When the vessel diameter can be varied but warning 
before fracture is given, in the form of a leakage i.e., 
the crack grows so that the flaw depth at failure is 
atlcnsb as much as the vessel wall thickness. In either 
of the two cases, the design should be optimal with 
respect to material cost, service life and the weight 

of the structure. The utility value of the design is 
considered as the net return from the design, viz, the 
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return from service life leas that incus***! -to- initial 
cost and the maintenance cost dependent on the weight of 
the structure. The material that maximizes the utility 
is the optimal material. The utility Up, the net 
monetary return is taken as given below: 

11 1 = V W ~ C ..(9.1) 

T Jp = net monetary return in Rupees 
ITg = service Tide in cycles 
w = weight of the vessel in lbs/inch length 
o = material cost Rs/inch length 

* constant denoting the return from one cycle of 
service = 0.5, in this example 

= instant denoting the cost incurred per cycle 

of service defending on weight = 0.01 in this example 

= denoting •••'» the total length of the vessel = 1000 in. 
in this example 

The return from service life and cost incurred as assumed 
in the o resent problem are shown in Figs. 9.2 and 9.3. 

It might be noted that in its totality, the problem of 
optimal design with material choice is a multifactor 
optimization problem (4, 123), wherein an optimum with 
reference to different criteria viz, weight, cost and 
service life is sought. It is necessary that a suitable 
trade off he made, so that the value of the design or 
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utility is maximized, The -definition of utility as is 
given in uhe present problem is not necessarily the most 
general one and any other alternate definition could as 
v/rll be given depending on the designer’s knowledge of 
the outcomes of various design decisions. An alternate 
form of utility model with trade off between various 
design criteria is attempted by Rake (5) and Murthy (6); 
however' further work is necessary to apply such models 
to practical design problems. In the present problem 
the utility as defined in Eqrt. 9*1 is used as a measure 
of optimality in the material choice. 

9 , 3 - late rials C ons idered and theirirPropcrties 

Six ” candidate” materials are considered (for 
deciding the best one) and their properties aro given 
in 'Table 9.1. The materials are designated by the 
corresoonding serial number in Table 9.1, in the subse- 
quent discussion. Because of the nonavailability of 
various properties for the same material, some of the 
properties have been assumed and are indicated by an 
asterisk in the table. The rate of propagation of a 
part— through crack which is of interest in estimating 
the service life does not appear to have been studied 
for any material. Tetelman and McEvily (128) reported 
the rate of crack propagation in steel specimens on 
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TABLE 9,1: PROPERTIES OP MATERIALS CONSIDERED FOR 
PRESSURE VESSEL DESIGN 


SI* 

No. 

Material 

Yield 

stre- 

ngth 

ksi 

Ulti- 

mate 

stre- 

ngth 

ksi 

Fract- 

ure 

tough- 

ness 

ksi 

fin 

Densi- 

. * 

1 7 ,, 

lh/cu 

Youngs 

modu- 

lus 

10 6 

psi 

Oost 

Rs. 

per 

lt>* 

Stra- Crack 
in at propa 
neck- gati- 
ing on 

instar> fact- 
hilitvy.o^ P 

1 

1 

I 8 M-C 0 - 
Mo steels 
(t?6) grade 
300-1 

286 

345* 

51.75 

0.283 

30 

3 

11 

1888 

2 

?50 

259 

345* 

68.4 

0.283 

30 

4 

11 

1865 

3 

300-2 

242 

* 

345 

84.5 

0.283 

30 

5 

11 

1850 

4 

Aluminium 
Alloys 
(125 ) 

2014-^651 

64.3 

69.7 

28.2 

0.107 

10,6 

2 

11.6 

128 

5 

2219-187 

57.6 

69.3 

27,9 

0.107 

10.6 

2 

10,5 

246 

6 

707 5-^6 

73.2 

80.7 

25.8 

0.107 

10,6 

2 

11.5 

157.5 


* Assumed 



7075-T6 Aluminium alloy and AM 355 CRT steel. The rate 
of crack propagation is expressed as a function of 


modified stress intensity factor denoted by 


m 


V 


_g q_'c 

°TT 


) 0 J) } i/4- 
> ~U U * i 


in 1/2 x 10 2 


(9*2) 


where 


T r 


i(— --vr 

* modified stress intensity factor 


li 

a = nominal stress applied 
G 

c = initial crack length 
a - yield stress 




u 


ultimate stress 

strain at necking instability 


E = Young* s modulus 

The relationship between K m and rate of crack propagation 
ns given by Tetelman end McEvily ' is reproduced in 
?ig. 9.4 and is assumed valid for the high strength 
steels and Aluminium alloys considered in Table 9.1 
in the presence ofa part-through crack. The quantity 
(3 for various materials denoted as crack propagation 
factor in Table 9.1 is the expression in the denominator 
of Eqn. 9.2 and is given by 


P = c ) 
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1/4 


.. (9.3) 
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Th.e computed values of 8 according to Eqn. 9*3 are given 
in the last column of Table 9,1, 


3*3.4 Design Procedure with Emphasis on Material Choice 
The stress analysis of the problem of part-through 
crack in a curved sheet is not yet solved. The corres- 
ponding solution for a flat sheet relating the flaw shape 
to the fracture stress with correction for plasticity at 
the crack tip is given by Tiffany and Masters (65). 

The fracture toughness K for the case of part-through 
crack is given by 


E 


c 



« » 


(9*4) 


where K c = fracture toughness 
o = applied stress 

a, = flaw depth at the onset of fracture 

G Jf 

Q = flaw shape parameter at the onset of 

G X 

fracture and is given by 


Q = * 2 - 0.212 ( X -) 2 .. (9.5) 

7 

where 4> = elliptic integral of the second kind 

tc/2 2 2 0 

= / /{I - ) Sin^ 9} de . (9.6) 

o a 

a = depth of the flaw 

o — dfltn'nia i or axis of the flaw 



Tlie flaw shape parameter curves relating Q and (a/2c) 
are plotted for various stress levels (65) and are 
reproduced in Fig, 9.5. Knowing the fracture toughness 
n t a material Eqn. 9.4 enables one to know the critical 
o rack size that leads to fracture at a given stress 
level. For the various materials considered in Table 
9 . 1 , the relationship between stress a and the critical 
oracle size calculated according to Eqn. 9.4 ar e plotted 
and are shown in Figs. 9.6 to 9.11. To proceed further 
with the .design, an estimate of the initial crack size 
is required. The most effective cracks in the cylindrical 
vessel are those which are located longitudinally on 
the shell so that they open out under hoop tension due 
to internal pressure. Consequently, it is assumed that 
the critical crack to be considered is under tension due 
to hoop stress. To estimate the initial crack size, 
proof loading upto 90% of the yield stress level is 
resorted to and if the vessel has not failed during the 
loading, one can infer that the initial crack size is 
less than the value a cr /Q cr corresponding to the 90% 
stress level on the curve relating stress and critical 
crack size for a material. As a conservative estimate, 
the value of a cr /Q cr corresponding to 90% of the yield 
stress can be taken as the initial crack size. Since, 
the working stress is to be less than 60-® of the 
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considered. Since the nominal working stress should not 
exceed 60“ of the yield stress the thickness ' t< of the 
vessel is fixed on the basis 


whore 


t 

n 

® y 

D 


= — 

2 x o.ey 

- thickness in inches 
= internal pressure = 2000 psi 
= yield stress 
= diameter = 50" 


.. (9.7) 


’res o" critical crack sizes corresponding to 90% 

>n*’ 60% of the yield stress level corresponding to various 
V’t. rifils are obtained from Pigs. 9.6 to 9.11. fable 9.2 
gj vox the details of working stress, proof stress, the 
thickness of the vessel wall according to Eqn. 9.7 and 
the vain os of a /Q ot> at working and proof stress levels. 
Oorresr ending to the ratio of a/2e as assumed i.e. , 0.4 
the values of are determined from the Pig. 9.5 giving 

Cr J? 

the detail s of flaw share parameter curves at various 

stress levels. Knowing the values of Q cr corresponding 

to proof and working stress levels, i.e. the initial and 

final values of a cr , from Pig. 9.5 the individual values 

nr u 2c , and a.*- , 2c,,. are determined 

of ^initial’ ^initial final’ final 

and are given in fable 9.3. The values of the modified 
stress intensity factor corresponding to 60% working stress 
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I 

IAHLE 9.2: PRESSURE! VESSEL DESIGN DETAILS (SAEE LIEE DESIGN) 


Mater- 

ial 

forking 
stress 
0.6 0 ^. 

ksi 

Proof 
stress 
0.9 a y 

ksi 

Vessel 
walls 
thick- 
ness’ t' 
inches 

a cr/^cr 

Proof 

stress 

level 

7/o r king 

stress 

level 

1 

171 

256 

0.292 

0.012 

0.026 

2 

155.5 

233 

0.321 

0.026 

0.052 

3 

145 

218 

0.344 

O 

• 

O 

0.08 

4 

38.6 

57.7 

1.295 

0.08 

0.16 

5 

34.6 

51.9 

1.445 

0.08 

0,26 

6 

43.9 

66.6 

1.14 

0.02 

0,10 


5 


/ 
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TABLE 9. 3: CRACK DIMENSIONS AND RATES OE PROPAGATION 


Mater- 

ial 

a. * 
ini- 
tial 
incli- 
ne 

a final 

inches 

2c. . 

ini- 

tial 

inches 

2o fi- 

nal 

inch- 

es 

Modified 

stress 

Intensity- 

factor 

Km 

Rate of 
crack 
propa- 
gation 

Service 

life 

cycles 

1 

0.0224 

0.0527 

0.056 

0.132 

1.5x10 2 

=r10“ 7 

76x10 4 

2 

0.0485 

0.1055 

0.121 

0.262 

2.12x10 2 

6.5x10“ 6 

2.7x10 4 

3 

0.0745 

0.162 

0.1865 0.405 

2.4 xIO 2 

4 x10“ 6 

5. 49x1 0 4 

4 

0*1495 

0.324 

0.374 

0.81 

12.9 xIO 2 

10 -1 

very low 

5 

0.1495 

0.527 

0.374 

1.29 

6.05x10 2 

1 . 5x10~ 4 

6. 3x10 4 

6 

0.0374 

0.2027 

0.0935 0.509 

0.6 xIO 2 

a 10“ 7 

415x10 4 
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level are determined, from Eqn. 9.2 and the corresponding 
rates of oracle propagation are obtained from Pig, 9.4, 

ilo values of K and the rate of crack propagation 

A 2o . 

A T *fcj ~ 1 ' re 1'atle 9.3. The numher of cycles 

r'-nuired to grow the crack from 2c. ... _ to 2c^.. .. is 

initial final 

determined and is given in the last column of Table 9.3. 
tho procedure is repeated for various materials and the 
Tinal details are given in Table 9.4, along with the wei- 
ghts and material costs computed using the properties 
given in Table 9.1. The utility ^calculated according 
to Eqn. 9.1, for the various materials' choice is also 
given, in the last column of Table 9.4. It can be seen 
from a comparison of the net returns from various design 
decisions, the materials one and six yield a "profit" , 
while the rest involve a "loss". The 4th material is not 
evaluated for its utility since the service life is too 
low. Obviously material six Is the best choice as it 
yields maximum return and satisfies all other design 
criteria. As could be seen from Table 9,. 3 the crack 
depths at fracture given by a final are less than the wall 
thicknesses in all the cases indicating that fracture 
occurs without any prior indication. For the optimal 
design chosen viz, using material six with wall thick- 
ness 1.14 inches, the vessel should be put out of service 
before 415 x 10 4 cycles unless Its further life can be 
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TABLE 9 *.4: DESIGN SOLUTIONS USING VARIOUS MATERI.'J/S 

(SAFE LIFE DESIGN) 


Mater- 
x al. 

Vessel 
Wall 
Thick- 
ness 't> 
inches 

Weight 
of the 
vessel 
w per 
inch 
length 
lhs 

Material 

cost 

C per 
inch 

1 ength 

Service 

life 

N 

cycles 

Net 

return 

d r 

Rs. 

1 

0*292 

13.0 

39 

76 x 10 4 

243000 

2 

0*321 

14.3 

57.2 

2.7 x 10 4 

-47450 

3 

0.344 

15.3 

76.5 

5.44x 10' h 

-57500 

* 

4 

1.295 

21.8 

43.6 

< 10 

- 

5 

1.495 

24.2 

48.4 

6.3 x 10 4 

-32000 

6 

1.14 

19.2 

38.4 

415 x 10 4 

1248600 


* Not considered in comparative evaluation since the 

service life is very low. 


ensured by suitable nondestructive tests and direct flaw 
observations. 


9.3.6 Ma terial Choice for Fail-Safe Design 

TTiile designing for safety at the onset of fracture 
keeping all the conditions of design as above, viz, the 
working and proof stress levels and the materials, the 
shell diameter cannot be a priori fixed and has to be 
fixed along with the thickness so that the flaw penetrates 
through the full depth before fracture occurs. The 
quantity U cr /Q cr ) initial as well as final, in Table 9,2 
remains valid even in this case except that the wall 
thickness should not be more than a cr so that, the crack 
fully propagates through the wall before the vessel 

fractures; i.e. 

= t for fracture safe design •• (9.8) 

' 'cr 

v/htfro = wall tTiictosss 

a = flaw depth at fracture 

cr 

Since the value of Q cr is mown as before, from the 
Pip. 9.5 giving the flaw shape parameter curves for 

stress -ex, t is computed for — • 

■The final values of (a cr /Q cr ) and the computed vessel 
thickness for varies materials are given in Table 95. 
Since the working stress level is fixed as 60% of the 
yield stress, the corresponding diameter h of the vessel 



TA3DE 9.5: DESIGN SOLUTIONS USING YAEKXJS MATERIALS 

(fail - s:^e design) 


Mater- 

ial 

for 

0 

or 

thick- 

ness 

t 

inches 

Diame- 
ter D 
inches 

Service 
life N 
cycles 

Weight 
W Its/ 
in 

1 ength 

Cost 

C Rs. 
per 
in 

1 ength 

Utility 

Net 

return 

Rs. 

1 

0,026 

0.053 

8.55 

76x10 4 

0.40 

1.2 

376800 

2 

0.052 

0.105 

16.3 

2.7x1'0 4 

1.53 

6.12 

6980 

3 

0,0:' 

0.162 

23.8 

5. 44x1 0 4 

3.42 

17. 10 

8300 

■* 

4 

0.16 

0.324 

12.5 

< 10 

- 

- 

- 

5 

0,26 

0.53 

18.3 

6’. 3x10 4 

3.23 

6.46 

22940 

6 

0.10 

0.202 

8.9 

415x10 4 

0.60 

1.2 

2038800 


* Not considered in comparative evaluation since the 
service life is very low 



is calculated, knowing the vessel thiclaiess from the 
eauation 

2. t. 0.6 cr 

D = — — — * .. (9.9) 

2000 

where 1 = diameter in inches 

t = thickness in inches 
o = the vi eld stress in psi 

The value of the diameter so calculated for various 
materials is also given in Table 9.5. Since the material 
stress levels and flaw shape are the same as in the 
earlier case of safe life design, the service lives as 
calculated previously remain unaltered even in the case 
of fail safe design. The weight and cost corresponding 
to the design solutions along with the associated utili- 
ties computed according to Eon. 9.1 are also given in 
Table 9.5. Comparing the utilities corresponding to 
various design solutions, once again it can be noted 
that the sixth material is the optimal as it yields 
maximum utility or return. 

9.4 DISCUSSION 

The method of optimal material choice for fracture 
resistant design discussed In this chapter enables one to 
choose a material out of an available set of materials 
so that the utility of design solution is maximized. 
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Even though the utility is considered synonymous with, 
net monetary return from service life, cost incurred 
due to weight and initial cost, any alternative definition 
could he used if necessary for the objective function 
wherein more design criteria like reliability and defor- 
mation could also be included. It is worth noting that 
the optimality criterion and design variables considered 
in the present problem are different from those of Heer 
and fang ( 127 ), wherein they have attempted to find an 
optimum proof load level and the probability of failure 
so that the weight of the structure is minimized 
subject to total cost and reliability constraints. 

Also they consider the design problem with a single 
material, in contrast to the present problem where 
the number of materials are considered and proof load 
level is a priori fixed. 

Even though some of the material properties assumed 
in the problems solved and the effects of curvature of 
the vessel on crack propagation are considered negligible 
because of lack of information, the method discussed 
remains valid and can be used if more accurate infor- 
mation regarding the material properties and stress 
analysis of part-through crack on a curved surface are 
available. For example, fatigue failure is not 
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considered in the design of the vessels; if data regard- 
ing the number of cycles the material can withstand 
at various stress levels is available, the same should 
also be considered in the service life estimation. 
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CHAPTER ten 

S2?lAR^t. CONCLUSIONS AND REG CffiSENDATI ONS 
FOR EURD'iiER RESEARC H ' 


10.1 SUMMARY A1TD CONCLUSIONS 

In this investigation, some selected problems 
in the area of statistical aspects of size effects on 
strength and fracture of materials (i.e. brittle 
fracture statistics) are discussed. The object and 
scope of the thesis is presented in Section 1.5. A 
detailed discussion as relevant to the various problems 
studied, is given at the end of each chapter. In the 
following a summary and conclusions based on the 
present investigation are given, while potential areas 
of further investigation are given in Section 10.2. 

A critical review of the various statistical 
strength theories as applied to various materials has 
indicated the inadequacy of the existing approaches 
in characterizing the size dependent scatter in 
materials. The possibility of an alternative approach 
viz, constructing a strength distribution function 
suitable to a given material, as given by its size- mean stre- 
ngth relation is investigated. Following the worlc 
of Tsai and Kolslsy (87), the aforementioned problem 
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is formulated using the weakest link concept and a 
distribution function is obtained corresponding to the 
size-mean strength relation of Wei bull in a reparameter- 
ized form. Essentially the method of approach currently 
suggested is converse in nature to the existing methods, 
since the strength distribution function is obtained 
corresponding to the size-mean strength relation of a 
material ratber than assuming a priori the distribution 
function. Even though the strength distribution function 
is obtained in clo sod form for the case studied in the 
thesis, treatment of more general forms of size-mean 
strength relations needs numerical techniques. 

Even though concrete is not a material that 
obeys the weakest link concept, the pattern of decrease 
in mean strength and scatter with increase in specimen 
size suggests an idealization of the material as that 
obeying weakest link concept, for the purpose of character- 
ization of scatter in strength. Using the test results 
of ICadlec ek and Spetla (43), an attempt is made to 
characterize the size dependent scatter in direct tensile 
strength of concrete, by choosing the material parameters 
in the size-mean strength relation by trial and error. 

It is noticed that the experimental and theoretical 
strength distribution functions are in agreement with 
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varying degrees corresponding to different specimen 
sizes, indicating that the same material parameters 
are inadequate to characterize the scatter with equal 
accuracy at various specimen volumes. 

Subsequently some possible practical applications 
of statistical strength theories are studied. Because of 
the existence of size dependent scatter in strength 
of materials, a problem of interest is to study the 
minimum number of specimens required to be tested (as 
affected by specimen size) so that the mean strength 
of the material is predicted with a given accuracy 

V 

und reliability. Using the test results of ICadlecek 
\/ 

and Spetla and applying Student's t distribution, an 
attempt is made to study the variation in the minimum 
number of specimens to be tested (for a given error of 
five percent with ninety percent probability, in 
predicting the mean direct tensile strength of concrete), 
with the specimen volume. As it is observed that 
smaller number of larger specimens are required, an 
attempt is made to find the optimum specimen size 
and sample number so that the total cost of testing 
is a minimum. After a study of some test results, , it 
is found that more number of tests are required to 
predict tensile strength of concrete, than the compressive 
strength. 
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A unified phenomenological theory of size effects 
on strength, stiffness and ductile-brittle transition 
in concrete-like materials is developed. Concrete-like 
materials are considered to he two phase composite 
materials , ''■nth one of the phases being voids due to 
microcracking. An estimate of the voids due to micro- 
crocking is made in terms of the stress and size 
dependent probability of failure. Using an approximate 
expression for the effective Young's modulus of a two 
phase composite material in the presence of voids, an 
expression for the size and stress dependent Young's 
modulus is derived. r l’he expression and the underlying 
theory arc shown to predict higher stiffness and greater 
ductility associated with specimens of smaller volume. 
The predictions are observed to be qualitatively in 
agreement with experimental observations in which steel 
fiber reinforced cement mortar briquettes of different 
sizes are tested in tension. 

As an application of fracture mechanics in 
engineering design, the design problem of a pressure 
vessel in the presence of a part-through crack with 
emphasis on material choice is considered. Six candi- 
date materials are considered. The utility of the 
design decision or material choice is considered to be 
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same as the net monetary return, viz, the return from 
service life, less the cost of material end the cost 
incurred due to its weight in service life. 

10.2 RECCMMSHUATIOUS FOR FURTHER RESEARCH 

The formulation and the method of solution of 
the problem of size effects on strength of materials 
discussed in the thesis, needs further theoretical 
studies regarding the uniqueness of distribution 
functions corresponding to a given size-mean strength 
relation. Also closed form solutions and numerical 
solutions to obtain the distribution functions relevant 
to various forms of size-mean strength relationships 
(as suitable to different materials) ar° to be studied 
in detail so that the present formulation can be applied 
for different materials for engineering use. 

Analytical and experimental studies on various 
materials for finding the size dependent coefficient of 
variation under general stress states would be of 
value in standardizing testing methods. Such studies 
in conjunction with the studies of actual mechanisms 
of failure would also enable expressing the scatter 
in the strength of materials in terns of the character- 
istic material parameters. 

The problem of the effective elastic moduli of 



212 


a composite material in the presence of voids, (in 
o articular interface discontinuities "between various 
phases) will "be helpful in studying materials like 
concrete in a more rational manner. Also experimental 
studies on such materials for size effects and volume 
changes will he valuable in furthering the knowledge 
of the mechanics of such materials. 

The problem of a part -through crack in a curved 
sheet for stress analysis as well as for rate of crack 
propagation would be of immense value in engineering 
design of pressure vessels etc. 
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